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1.0  PROGRESS  ANO  TECHNICAL  REPORT  SUMMARY 

The  objective  of  this  program  is  to  perform  fundamental  studies  on 
SAW  resonators  leading  to  improved  bias  stability  characteristics  for  acceler¬ 
ation  sensors.  The  end  goal  is  to  develop  small,  low  cost  and  moderately 
accurate  accelerometers  for  tactical  guidance  systems. 

During  the  first  year  progress  has  been  made  in  several  areas  which 
are  discussed  in  this  report.  A  method  of  fabricating  dual -resonator  crystals 
with  absolute  and  differential  aging  characteristics  has  been  developed. 

Using  a  unique  method  of  measuring  process  dependent  shifts  in  resonator 
frequency  with  time,  both  absolute  and  differential  bias  stability  has  been 
studied.  In  conjunction  with  a  study  of  crystal  processing  methods  a  study  of 
the  surface  chemistry  of  electroded  quartz  resonators  using  such  techniques  as 
Auger,  ESCA,  and  mass  spectroscopy  was  carried  out.  The  quantitative  results 
of  these  studies  have  been  combined  with  frequency  stability  measurement  data 
to  arrive  at  a  mathematical  model  which  explains  long  term  frequency  shifts  in 
terms  of  logarithmic  rate  laws  associated  with  chemisorption,  oxidation  and 
stress  relief  in  surface  wave  resonators. 

Hybrid  RF  circuitry  for  prototype  dual  resonator  accelerometers  has 
been  developed.  Preliminary  results  using  dual  resonator  crystals  show 
extremely  good  stability,  typically  less  than  10"^  for  1  second  averaging 
times.  Noise  levels  due  to  electronic  circuitry  are  consistent  with  a  SAW 
sensor  dynamic  range  of  1x10®.  Tests  were  also  performed  on  integrated  noise 
levels  for  simulated  guidance  system  mission  times.  For  a  mission  time  of  3 
min  the  cumulative  displacement  error  rate  was  typically  less  than  1  meter/hour 
and  for  simulated  mission  times  of  20  min  the  cumulative  displacement  error 
rate  was  less  than  20  meters/hour.  This  data  is  to  be  compared  to  a  typical 
error  rate  for  a  moderately  accurate  accelerometer  of  1  nautical  mile  per  hour 
or  1851  meters/hour.  These  results  are  encouraging  and  during  the  second  year 
further  studies  of  SAW  accelerometers  will  be  performed. 
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2.0  TECHNICAL  PROBLEM 

Accelerometers  are  transducers  which  convert  sensed  acceleration  into 
analog  or  digital  signals.  They  are  used  in  a  wide  variety  of  applications, 
such  as  in  sensing  seismic  and  vibrational  motions;  but  the  most  stringent 
requirements  for  these  devices  occur  in  their  applications  for  guided  missiles. 

The  development  of  guided  missiles  in  the  1940* s  required  the 
development  of  accurate  accelerometers  for  steering  and  stability.  This  led 
to  the  development  of  accelerometers  for  autonomous  inertial  navigation  of 
missiles,  aircraft  and  ships.  They  are  also  used  for  autopilots  and  smoothing 
for  navigational  mechanizations  using  satellite  data. 

The  technological  approaches  for  accelerometers  of  the  1940's  are 
still  in  use  today.  These  include  integrating  gyroscopic  accelerometers  of 
the  type  used  on  the  German  V- 2  rockets  (R-l).  Refinements  to  these  approaches 
have  improved  these  instruments  in  many  respects,  but  they  have  not  made  them 
less  expensive.  They  require  precision  machining  and  assembly  of  many 
mechanical  parts,  which  cannot  be  accomplished  at  low  cost.  Typical  unit 
costs  for  accelerometers  used  in  inertial  navigation  range  from  a  few  thousand 
dollars  to  tens  of  thousands  of  dollars. 

There  is  a  constant  need  for  less  costly  accelerometers,  especially 
for  applications  such  as  cruise  missiles  and  tactical  missiles,  for  which 
large  numbers  of  accelerometers  are  required.  Many  technological  approaches 
have  been  proposed  for  meeting  this  need.  These  include  vibrating  string 
(R-2),  laser  (R-3),  piezoresistive  (R-4),  piezoelectric  (R-5),  and  acoustic 
(R-6,7)  accelerometers.  All  of  these  approaches  use  a  'proof  mass'  which  is 
connected  to  the  host  vehicle  by  means  of  a  stress  sensor.  The  stress  sensor 
measures  the  force  applied  to  the  proof  mass,  which  is  proportional  to  the 
sensed  acceleration  of  the  host  vehicle.  These  are  'open  loop'  sensors, 
because  the  force  is  not  applied  by  a  precision  servo  loop,  as  in  the  case  of 
the  gyroscopic  and  pendulous  accelerometers. 
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The  surface  acoustic  wave  (SAW)  resonator  transforms  strain  input  to 
frequency  output  (R-8,10).  When  used  as  an  open-loop  accelerometer  without  a 
rebalance  servo  loop,  the  SAW  oscillator  functions  as  a  solid-state  digital 
accelerometer  with  the  desirable  characteristics  of  low  power  consumption,  low 
weight,  small  size,  and  potentially  very  low  cost.  In  order  to  meet  perform¬ 
ance  requirements  for  moderately  accurate  inertial  navigation  systems,  the 
frequency  stability  must  be  high.  Because  SAW  resonators  operate  well  into 
the  UHF  range,  SAW  sensors  are  normally  configured  as  dual  resonator  controlled 
oscillators  and  one  resonator  acts  as  a  frequency  reference  or  standard.  In 
this  case  the  differential  frequency  stability  determines  the  noise  level  of 
the  accelerometer  output. 

The  objective  of  this  program  is  to  perform  fundamental  studies  on 
SAW  resonator  processing  methods  that  yield  stability  characteristics  which 
successfully  meet  the  bias  stability  requirements  of  inertial  navigation 
accelerometers.  The  technical  problems  involved  are  SAW  resonator  fabrication 
methodology,  surface  chemistry  of  electroded  single  crystal  quartz,  absolute 
and  differential  aging  measurements,  design  and  characterization  of  high 
frequency  oscillator  circuitry,  and  common  mode  noise  rejection  including 
temperature  effects.  During  the  first  phase  mathematical  models  for  predict¬ 
ing  and  explaining  frequency  shifts  in  SAW  resonators  configured  as  acceler¬ 
ation  sensors  were  compared  with  experimental  data.  Statistically  valid 
measurement  methods  were  developed  in  conjunction  with  improved  SAW  sensor 
crystal  fabrication  methods. 
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3.0  TECHNICAL  RESULTS 

During  the  first  year  of  this  program,  the  specific  technical  tasks 
were  the  following: 

1.  To  develop  methods  for  experimentally  evaluating  single  and 
differential  bias  stability  in  two-resonator  crystals.  This  problem  involved 
the  use  of  automatic  data-acqui sition  hardware  and  software  for  measuring  long 
term  frequency  shifts. 

2.  To  develop  mathematical  models  for  predicting  frequency  shifts  in  SAW 
resonators  based  upon  the  known  behavior  of  physical  processes.  This  model 
was  used  to  interpret  the  experimentally  measured  aging  data. 

3.  To  perform  studies  of  the  chemical  interactions  which  take  place 
within  and  on  surface  wave  structures  when  sealed  under  specified  environ¬ 
mental  conditions.  The  problem  is  to  explain  the  kinematics  of  chemical 
interactions  in  sealed  crystal  packages  and  relate  these  effects  to  the 
surface  chemistry  of  the  quartz  crystal  and  metallic  surface  wave  electrode 
structure. 

4.  To  design  and  construct  hybrid  rf  circuitry  for  prototype  SAW 
accelerometers.  The  circuitry  must  be  low  in  cost  and  power,  and  provide  long 
term  stability.  The  technical  problem  was  to  operate  two  oscillators  from  a 
single  crystal,  without  experiencing  lock-up  of  the  oscillators  to  a  common 
frequency. 

5.  To  assess  the  inertial  errors  of  prototype  SAW  acceleration  sensors 
in  terms  of  the  integrated  velocity  and  displacement  errors.  This  required 
integration  of  the  SAW  sensor  digital  output  twice,  once  to  derive  velocity 
and  once  to  derive  displacement  from  velocity.  The  inertial  errors  were 
studied  as  a  function  of  mission  time  at  constant  temperature. 
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3.1  SAW  Resonator  Fabrication 

The  surface  acoustic  wave  resonator  is  a  planar  electrode  structure 
photol ithographical ly  defined  on  a  suitable  piezoelectric  material.  The  SAW 
resonator  functions  in  much  the  same  manner  as  a  conventional  quartz  oscil¬ 
lator  crystal.  Acoustic  waves  are  confined  within  a  cavity  whose  boundaries 
are  accurately  maintained.  The  Q  of  the  cavity  is  determined  by  the  material 
losses  and  cavity  leakage.  The  same  excellent  frequency  control  properties  of 
bulk  wave  quartz  resonators  is  also  achieved  by  the  SAW  resonators  provided 
the  fabrication  techniques  are  closely  controlled. 

The  planar  electrode  pattern  of  a  SAW  resonator  contains  an  inter- 
digital  transducer  which  provides  electrical -to-mechanical  as  well  as 
mechanical -to-electrical  signal  transduction.  Shown  in  Fig.  1  is  a  photograph 
of  the  resonator  geometry  used  in  this  study.  The  interdigital  transducer 
contained  40  electrode  finger-pairs  with  a  periodicity  of  8.3  microns.  On 
either  side  of  the  input/output  transducer  are  two  gratings  which  contain  400 
reflecting  stripes  with  a  periodicity  slightly  less  than  in  the  transducer. 

The  gratings  act  as  mirrors  when  the  acoustic  wavelength  is  approximately 
equal  to  twice  the  grating  periodicity.  In  this  frequency  range  all  of  the 
surface  acoustic  wave  energy  is  confined  within  the  cavity  formed  by  the  two 
gratings. 

The  complete  SAW  resonator  fabrication  process  used  in  this  study  is 
shown  in  Fig.  2.  The  starting  material  is  Y-cut  single  crystal  quartz  purchased 
from  Valpey-Fi sher  with  an  ultrafine  polish  on  one  surface.  Typically  the  rms 
surface  roughness  is  less  than  50  angstroms.  The  first  row  of  Fig.  2  represents 
the  photolithographic  process  used  to  define  the  SAW  electrode  pattern  of  Fig.  1. 

A  photomask  with  over  100  resonators  covering  a  two  inch  square  area  is  used  to 
define  and  etch  the  electrode  pattern  in  a  vacuum  deposited  aluminum  film.  In  our 
work  three  types  of  metalization  systems  were  used,  (1)  electron-beam  evaporation 
of  chrome  (50  A)  and  aluminum  (1000  A),  (2)  sputter  deposited  pure  aluminum 
(1000  A),  and  (3)  sputter  deposited  copper-doped  (4%)  aluminum  (1000  A). 
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Fig.  1  Photomicrographs  of  SAW  resonator  structure  used  in  absolute 
and  differential  aging  study. 
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After  the  metal  electrode  pattern  has  been  delineated  the  metal  thickness  is 
evaluated  for  the  desired  value  as  well  as  a  check  of  metal  thickness 
uniformity.  Uniformity  of  the  deposited  metal  was  important  to  achieving 
resonators  closely  matched  in  frequency.  An  rf  probing  technique  was  used  to 
measure  the  resonator  frequencies  prior  to  any  other  processing  of  the 
wafers.  The  frequency  sensitivity  of  aluminum  electrode  thickness  was 
typically  1.67  kHz/A  at  nominally  375  MHz.  In  the  early  fabrication  runs  the 
frequency  variation  was  typically  1000  ppm.  It  was  necessary  to  the  maintain 
a  deposition  uniformity  of  less  than  5%  in  order  to  achieve  runs  with  less 
than  200  ppm  frequency  variation. 

In  order  to  achieve  high  Q  in  SAW  resonators  efficient  surface  wave 
reflection  in  the  grating  regions  of  the  resonator  must  occur.  This  is 
accomplished  in  the  next  step  (second  line  in  Fig.  2).  A  photoresist  (AZ1350) 
is  spun  onto  the  wafers  and  patterned  such  that  windows  are  opened  in  the 
reflecting  grating  regions  of  the  resonator.  Following  this  a  freon  (CF4) 
reactive  ion  etch  of  the  quartz  substrate  material  was  performed.  The  photo¬ 
resist  prevented  any  etching  of  the  transducer  regions  and  the  aluminum  in  the 
grating  regions  acted  as  an  etch  mask  allowing  only  the  quartz  between  the 
gratings  to  be  etched.  This  produced  grooves  between  these  electrodes  which 
enhanced  the  surface  wave  reflection  at  the  groove  edges  and  also  lowered  the 
resonant  frequency  of  the  gratings  to  coincide  with  the  transducer  frequency. 
This  dispersion  was  accounted  for  by  designing  the  photomask  with  different 
periodicity  in  the  gratings  and  transducer.  After  plasma  etching  the  resist 
was  removed  and  the  resonators  again  rf  probed.  A  rework  cycle  was  used  to 
place  the  resonators  within  a  prescribed  frequency  range.  Typically  the  total 
frequency  change  due  to  plasma  etching  was  1%  and  resonator  Q's  were  nominally 
10-15,000  at  375  MHz. 

After  plasma  etching  the  wafers  were  diced  and  probed.  The  dicing 
procedure  was  designed  to  yield  chips  with  two  resonators  on  the  same  chip. 
This  structure  was  termed  the  'dual  crystal'  resonator  and  was  used  exten¬ 
sively  in  our  testing  of  stability  and  aging.  After  dicing  the  chips  were 
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mounted  using  a  polyimide  epoxy  which  was  cured  at  350°C.  Following  mounting, 
typically  to  a  TO-5  header,  the  individual  resonators  were  wire  bonded  and  a 
pre-test  of  the  resonator  characteristics  was  performed.  Pre-testing  of  the 
resonators  yielded  device  frequency,  Q,  and  resistance  before  sealing. 

Sealing  was  performed  in  a  controlled  glove  box  environment.  The 
controls  included  an  automatic  pressure  controller  and  a  gaseous  dri-train 
system  for  removal  of  all  water  vapor  and  oxygen  from  the  glove  box  enclosure. 
A  capacitance-discharge  resistance  weld  was  used  to  seal  the  TO-5  crystal 
enclosures.  Input  and  output  vacuum  load-looks  were  used  to  transport  the 
resonators  through  the  sealing  system.  The  input  load-look  chamber  was  also  a 
vacuum  oven  and  a  24  hour  pre-seal  bake  at  200°C  was  routinely  performed. 

After  this  vacuum  bake  but  prior  to  actual  sealing,  each  crystal  was 
calibrated  to  a  desired  frequency  using  a  short  plasma  etch  in  freon.  The 
calibration  system  was  installed  inside  the  glove  box  and  is  shown  in 
Fig.  3.  The  etching  geometry  is  of  the  parallel -plate  type  with  the  SAW 
resonator  attached  to  the  anode  and  rf  power  applied  to  the  cathode. 

Provisions  for  monitoring  the  resonator  frequency  were  made  using  rf  feed¬ 
throughs  to  the  vacuum  chamber.  The  SAW  resonators  could  be  tuned  by  as  much 
as  1000  ppm  without  degradation  and  could  be  placed  to  within  +10  ppm  of  the 
desired  frequency. 

The  sealing  operation  itself  was  carried  out  under  a  vacuum  of  10"7 
torr.  After  sealing  the  devices  were  bombed  with  helium  and  leak  checked  with 
a  mass  spectrometer.  Those  devices  which  showed  no  leakage  were  then  put 
through  a  series  of  tests  to  determine  their  long  term  frequency  stability 
under  various  environmental  conditions. 

3.2  Surface  Acoustic  Wave  Resonator  Aging  Measurements 

The  quantitative  determination  of  resonator  aging  was  a  major  task  of 
this  study.  The  amount  by  which  a  quartz  resonator,  bulk  or  SAW,  changes  with 
time  is  very  small.  Measurement  accuracies  of  0.1  ppm  or  less  are  required 
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over  time  periods  ranging  up  to  one  year.  Because  the  resonator  is  a  passive 
device  the  conventional  method  is  to  build  oscillator  circuitry  with  the 
resonator  as  the  frequency  feedback  element.  However,  when  considering  the 
large  numbers  of  resonators  to  be  tested  this  approach  was  not  practical. 
Instead,  a  method  of  determining  the  individual  resonator  series  resonant 
frequency  periodically  with  an  automated  network  analyzer  system  was  developed. 
In  this  technique  absolute  frequency  accuracy  was  obtained  by  phase  locking 
the  network  analyzer  to  an  external  frequency  standard,  such  as  a  frequency 
synthesizer  locked  to  a  rubidium  clock.  In  practice  the  rubidium  standard  was 
not  needed  for  time  periods  less  than  1  year. 

The  system  which  was  first  built,  successfully  monitored  the 
frequency  stability  of  12  resonators  in  two  ovens.  Each  oven  was  fitted  with 
a  single  pole,  6-throw,  coaxial  switch  which  selected  which  resonator  was  to 
be  tested.  The  coaxial  switchs  were  driven  by  an  HP3495A  scanner  digitally 
programmed  by  an  HP9825A  desktop  calculator.  The  coaxial  switch  outputs  were 
connected  to  the  A  and  B  ports  of  an  HP  8505A,  0.1-1300  MHz,  network  analyzer, 
also  controlled  by  the  desktop  calculator  via  an  IEEE-488  data  bus.  The 
network  analyzer  was  phase  locked  to  an  HP8660C  frequency  synthesizer  which 
acted  as  the  system  frequency  standard  as  well  as  the  primary  scanning 
frequency  source  for  the  network  analyzer.  The  crystal  ovens  were  controlled 
at  a  temperature  of  85°C,  ±0.1°C.  All  of  the  aging  data  taken  in  the  first 
two  quarters  of  this  study  was  obtained  using  this  prototype  measurement 
system. 

Software  for  the  aging  measurements  was  designed  to  control  the  system 
hardware  and  take  measurements  periodically.  This  was  accomplished  by  using  a 
real  time  clock  within  the  desktop  computer.  Resonator  aging  data  was  not 
continuously  taken,  but  rather  a  system  of  interrupts  was  used.  Whenever  the 
real  time  clock  generated  an  interrupt,  the  controller  automatically  directed 
the  scanner  to  step  through  the  crystals  and  the  network  analyzer  to  determine 
the  scattering  parameters  of  each  crystal.  Impedance  characteristics  and  series 
resonant  frequencies  were  then  derived  in  software  and  stored  on  magnetic  tape. 
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Extensive  testing  showed  that  this  method  of  measuring  resonator 
aging  was  satisfactory  and  a  larger  oven  with  an  integral  rf  switching  matrix 
was  designed  for  monitoring  the  aging  of  36  crystals  simultaneously.  By  using 
interchangeable  PC  boards,  the  aging  of  either  18  TO-5  packages  with  resonator 
pairs,  or  36  packages  containing  single  resonators,  could  be  monitored. 
Resonator  groups  of  six  were  connected  to  a  single-pole  six-throw  (SP6T) 
coaxial  switch.  The  six  switches  were  connected  to  a  seventh  in  a  switching 
matrix  arrangement.  Figure  4  shows  a  photograph  of  the  oven  containing  the  PC 
board  into  which  the  TO-5  packages  were  inserted  and  the  coaxial  switching 
matrix. 

The  resonators  were  periodically  interrogated  through  the  switches  by 
a  network  analyzer,  frequency  synthesizer  and  counter,  all  of  which  were 
controlled  by  a  desktop  computer  as  shown  in  Fig.  5.  The  resonant  frequencies 
were  sequentially  measured,  the  oven  temperature  recorded,  and  the  time  logged 
from  the  computer's  real  time  clock.  In  this  manner,  the  frequency/temperature 
drift  was  tracked  with  time.  Software  computed  the  differential  drift  between 
packaged  pairs  of  SAW  resonators. 

Experimental  data  taken  with  the  36  crystal  switching  oven  is  shown 
in  Fig.  6.  The  plotting  of  experimental  data  was  controlled  by  interactive 
software  which  could  be  directed  to  simultaneously  plot  the  frequency  changed 
(ppm)  of  all  36  crystals  as  a  function  of  linear  time.  Fig.  6a,  or  as 
logarithmic  time  (log[T+l]),  Fig.  6b.  Examination  of  such  data  clearly 
indicates  the  anomalous  behavior  of  some  crystals  and  more  importantly  shows 
the  overall  or  average  aging  characteristics  of  the  entire  group.  Knowledge 
of  the  latter  was  invaluable  in  obtaining  a  better  understanding  of  the 
functional  dependence  fabrication  methods  play  in  SAW  resonator  aging. 

Breaks  in  the  aging  data  curves  represent  times  when  the  measurement 
system  was  interrupted  for  other  types  of  measurements  or  tests.  Overall 
stability  and  repeatibility  was  excellent.  Even  through  the  interrupt  time 
was  several  days  in  some  cases,  the  curves  continue  without  any  apparent 
discrepancies  in  aging  trends.  Individual  or  paired  plots  of  crystal  aging  or 
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6.  Aging  data  representative  of  one  group  of  36  resonator  crystals.  In 
Fig.  6a  the  absolute  aging  Is  plotted  vs  linear  time  and  In  Fig.  6b 
the  same  data  is  shown  vs  logarithmic  time. 
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scattering  parameters  could  easily  be  obtained  for  more  detailed  study. 
Examples  of  these  types  of  plots  as  well  as  measurement  system  software  are 
discussed  in  Appendices  I -IV  of  this  report. 

3. 3  Surface  Chemistry  Related  to  SAM  Resonator  Aging 

A  series  of  ultra-high  vacuum,  surface  analysis  experiments  were 
carried  out  on  the  bulk  materials  of  the  SAW  oscillators  to  evaluate  the 
effects  of  the  processing  steps  upon  the  aging  characteristics  of  the  devices. 
In  one  series  of  experiments  the  individual  processing  steps  were  investigated 
to  establish  the  level  of  impurities  which  they  contribute  in  resonator 
materials.  In  the  second  series  of  experiments  the  device  materials  were 
depth  profiled  by  secondary  ion  mass  spectroscopy  (SIMS)  to  establish  the 
extent  of  diffusion  of  the  aluminum  and  chromium  reflector  material  layers 
into  the  quartz  substrate  and  each  other. 

The  investigation  of  the  contamination  levels  associated  with  the 
bulk  materials  and  the  processing  steps  were  carried  out  using  the  techniques 
of  x-ray  photoelectron  spectroscopy  (XPS)  and  SIMS.  The  XPS  analysis  of  the 
surface  of  the  quartz  substrate  material  generally  showed  little  chemical 
reactivity  of  the  quartz  with  atmospheric  gases;  however,  the  SIMS  analysis 
identified  a  significant  involvement  of  water  in  the  top  20A  of  the  surface 
with  the  SiOH+  being  generated  as  a  secondary  ion  under  argon  ion  bombard¬ 
ment.  In  addition  to  the  presence  of  hydroxide  groups  on  the  quartz  surface, 
hydrocarbons  at  the  level  of  0.5  atom  %  were  also  found  to  be  present  at  the 
quartz  metal  interface. 

Samples  of  the  substrate  material  upon  which  a  thin  layer  (40  A)  of 
chromium  followed  by  aluminum  had  been  deposited,  were  also  examined  for 
incorporation  of  contaminants.  XPS  and  SIMS  analysis  showed  that  the  layers 
were  free  from  any  major  contaminants  with  only  potassium  and  water  being 
observed  at  a  level  of  about  0.1  atom  %.  These  levels  of  contamination  were 
compared  against  those  observed  in  single  crystal  aluminum  samples  exposed  to 
air  and  were  found  to  be  the  same  magnitude. 
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Depth  profiling  of  the  deposited  layers  of  Cr-Al  upon  the  quartz- 
substrate  substrate  were  carried  out  again  using  XPS  as  well  as  SIMS.  The 

profiling  studies  showed  the  layered  structure  to  be  composed  of  six  distinct 

regions  as  shown  in  Fig.  7.  The  upper  layer  of  aluminum  showed  no  major 
impurities;  however,  the  chromium  substrate  layer  had  diffused  into  aluminum 
to  a  distance  of  approximately  100  to  300  A.  The  chromium  layer  itself  was 
not  pure  metal,  containing  a  significant  fraction  of  oxide.  At  the  quartz- 
chromium  interface  the  SIMS  analysis  clearly  identified  the  presence  of 
hydroxides,  both  the  presence  of  the  SiOH+  ion  and  a  marked  increase  in  the 

amplitude  of  the  Cr+  ion  intensity,  at  the  interface  which  is  normally 

associated  with  the  presence  of  oxygen.  The  impurity  layer,  however,  was 
again  clearly  defined  and  quite  thin,  covering  only  about  10  to  20  A  at  the 
interface.  The  chromium  was  found  to  have  diffused  into  the  quartz  substrate 
to  a  distance  of  300  to  600  A,  showing  an  exponential  decrease  in  concentra¬ 
tion  with  sputtering  time  into  the  quartz. 

This  marked  diffusion  of  chromium  into  the  A1  layer  and  the  quartz 
lead  to  the  investigation  of  an  improved  reflector  material  made  of  layers  of 
A1  containing  about  4%  Cu.  This  material  was  investigated  on  samples  which 
had  been  specifically  designed  for  SIMS  analysis  and  which  had  been  exposed  to 
all  the  device  processing  steps.  The  quartz  substrates  were  cleaned  and 
metalized  with  the  Al/Cu  material.  Then  a  photoresist  was  laid  down,  developed, 
and  acid-etched  to  leave  a  pattern  of  Al/Cu  rectangles  on  the  exposed  quartz 
substrate.  The  samples  were  then  acetone-rinsed,  A-40  stripped,  and  nitric 
acid  cleaned.  Some  of  the  samples  were  then  exposed  to  plasma  etching  in 
freon  gas  (CF4  +  4%02). 

SIMS  analysis  of  the  metalized  areas  clearly  showed  the  presence  of 
the  Cu  impurity  in  the  A1 ;  however,  no  evidence  was  found  for  any  residual  Cu 
in  the  exposed  quartz  after  the  etching  and  cleaning  process.  The  A1  surface 
before  the  plasma  etch  showed  a  clear  A10x  oxide  film  on  the  surface.  This 
oxide  layer  was  quite  thin  (10-20  A)  and  the  bulk  elemental  A1  could  be 
clearly  distinguished  under  the  oxide  by  XPS,  which  observes  only  the  top  30 
to  50  A  of  the  surface. 
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Exposure  of  the  device  surface  to  the  CF^  +  O2  plasma  etch  clearly 
formed  a  surface  layer  of  CFX.  The  layer  was  approximately  20  to  30  A  thick 
and  had  a  fl uori ne-to-carbon  ratio  of  3  to  3.4.  A  marked  shift  in  the  A1  XPS 
lines  toward  higher  binding  energies  indicated  CFX  was  associating  with  the 
A1 ,  possibly  as  A10XCF3.  Samples  of  the  plasma  treated  devices  were  heated  to 
350°C  in  N2  and  showed  no  marked  decrease  in  the  intensity  of  the  XPS  signals 
from  the  fluorocarbon  surface  layer,  demonstrating  the  layer  to  be  chemically 
bound  and  quite  stable  rather  than  simply  adsorbed  gas.  Heating  the  samples 
in  air  or  oxygen  as  an  annealing  process,  however,  removed  80  to  90%  of  the 
fluorocarbon;  returning  the  surface  to  a  dominantly  41 0X  layer.  The  residual 
oxide  layer  still  contained  a  small  quantity  of  the  A10XCF3  type  structure  and 
was  overall  slightly  thicker  than  the  original  oxide  film  prior  to  the  plasma 
etching  and  annealing  procedure. 

Studies  of  the  surface  chemistry  of  electroded  quartz  resonators 
using  techniques  such  as  SIMS,  Auger,  and  mass  spectroscopy  have  clearly 
indicated  that  chemisorption  and  oxidation  are  taking  place  and  that  inter¬ 
facial  and/diffusion  layers  do  form  in  these  electrodes.  This  information  was 
supported  by  actual  SAW  resonator  isothermal  aging  data.  In  the  case  of 
chrome,  the  aging  with  time  was  2-5  times  the  rate  with  pure  aluminum  elec¬ 
trodes.  The  aging  rate  was  even  lower  when  copper  doped  aluminum  was  used. 

In  addition,  degradation  under  high  rf  power  levels  was  less,  indicating  mass 
transport  was  less.  These  studies  have  contributed  greatly  to  our  understand¬ 
ing  of  SAW  resonator  long  term  aging. 

3.4  Absolute  and  Differential  Aging  in  SAW  Resonators 

An  important  objective  of  this  study  was  to  develop  mathematical 
models  for  predicting  frequency  shifts  in  SAW  resonators  based  upon  the  known 
behavior  of  physical  processes.  Three  types  of  aging  models  were  postulated 
for  these  processes, 
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f-  =  A  log  (BT  +  1) 

^  =  A(1  -  e’BT) 
r-=  A{T}8 

The  first  model  is  applicable  to  chemisorption,  oxidation,  and  stress-relief 
( R-l 1 ,R-1 2) .  The  second  applies  to  first  order  chemical  reactions  or  adsorp¬ 
tion/desorption  of  physisorbed  monolayers  upon  the  substrate.  The  final  model 
has  been  identified  with  diffusion  processes  when  B  is  near  1/2. 

In  virtually  every  experimental  SAW  resonator  which  had  been  sealed 
under  vacuum  without  significant  leakage,  the  aging  data  was  always  of  the 
form  predicted  by  the  first  model.  This  result  strongly  suggested  that  the 
aging  of  SAW  resonators  was  dominated  by  either  chemisorption,  oxidation, 
stress-relief,  or  a  mixture  of  all  three.  Physically  these  types  of  processes 
are  expected  and  were  clearly  detected  during  our  work  on  electroded  quartz 
surface  chemistry.  Oxidation  of  the  aluminum  electrodes  was  unavoidable. 
Stress  induced  by  crystal  mounting  as  well  as  electrode  formation  during 
evaporation  or  sputtering  was  also  present.  Chemisorption  was  not  easily 
detected,  however  diffusion  of  metal  ions,  particularly  when  chromium  was 
used,  was  confirmed  by  XPS  and  SIMS  spectroscopy. 

Theoretically  plotted  logarithmic  frequency  changes  based  on  the 
first  model  are  shown  in  Fig.  8a, b.  A  family  of  four  curves  are  shown  plotted 
vs  linear  time  in  Fig.  8a.  The  same  four  curves  are  plotted  vs  logarithmic 
time  in  Fig.  8b.  The  logarithmic  time  plots  show  that  after  a  short  initial 
time  the  curve  becomes  a  straight  line.  This  behavior  is  typical  of  the 
logarithmic  rate  law  and  is  closely  supported  by  data  (see  Fig.  6).  As  an 
example  consider  the  typical  device  aging  data  plotted  in  Fig.  9a.  In  Fig.  9b 
the  same  data  is  plotted  versus  four  theoretical  values  of  A  in  the  rate 
law.  The  best  fit  appears  to  be  A=4 . 25  ppm  and  B=0. 12/hr.  The  aging  rate  at 
the  end  of  one  year  for  this  crystal  can  be  calculated  from  the  derivative  of 
the  rate  law. 
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Fig.  9.  (a)  Typical  absolute  aging  data  vs  linear  time  for  a  SAW  resonator. 

The  same  data  Is  plotted  (b)  vs  logarithmic  time  and  compared  to  rate- 
law  models  with  different  A  parameters. 
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n  _  AB  log(e) 

R  “  "TTT+l 

Substituting  for  A  and  B,  the  aging  rate  at  the  end  of  one  year  is  0.000211  ppm/hour 
or  1.844  ppm/year.  Alternatively  the  time  required  to  achieve  a  given  aging 
rate  can  be  expressed  as, 

T  =  (AB  -  1)/B 

where  R  is  the  rate  and  T  is  the  required  time.  Substituting  the  appropriate 
data  for  the  experimental  crystal  shown,  the  time  required  to  reach  an  aging 
rate  of  1  ppm/year  is  1.8448  years. 

For  large  times  the  logarithmic  law  can  be  approximated  by 
£■  -  A  log(T) 

where  the  controlling  parameter  is  A.  The  rate  of  aging  using  this  approx¬ 
imation  is  simply, 

R  =  A  1 og(e)/T 

These  simple  relations  can  be  used  to  determine  the  value  of  A  needed  to 
achieve  a  given  degree  of  aging  after  a  prescribed  time.  For  example  to 
achieve  a  1  ppm/year  aging  rate  after  30  days  pre-aging  will  require 

A  =  (12  log(e))-! 

Assuming  the  constant  B  does  not  change,  this  curve  is  plotted  in  Fig.  9b. 

Also  shown  are  pre-aging  curves  under  the  same  conditions  for  60  and  90  days 
pre-aging.  Since  a  pre-aging  time  of  30  days  or  less  is  desired,  it  is  clear 
that  a  total  aging  after  30  days  of  less  than  0.5  ppm  is  required. 

Having  established  the  particular  rate  law  governing  SAW  resonator 
crystal  aging  a  study  of  differential  aging  was  made.  Differential  aging  in 
the  present  study  was  of  interest  since  it  directly  contributed  to  differential 
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bias  stability  in  SAW  sensors  using  dual  crystals  to  detect  acceleration.  The 
objective  was  to  predict  under  what  conditions  the  differential  aging  was 
least.  Assuming  two  resonators  with  the  aging  behavior  given  by. 


F1  •  Fl0d  +  A1  log(BiT  +  1)) 
F2  =  f20^  +  A2  1og(B2T  +1)) 


the  change  in  their  difference  frequency  is  given  by 


where 


Fi  nAl 

(B  T  +  1)  10  1 
AF  -  AF  =  log  - p— r- 

(B2T  .  1)  20  2 


AFo  =  F10  "  F20 


Two  conditions  are  likely  to  occur,  namely  A^  #  A2,  and  B1  #  B2.  These  are 
plotted  in  Figs.  10  and  11  for  hypothetical  crystals.  When  the  B  constants 
are  not  equal  a  small  change  occurs  initially  but  after  a  time  large  compared 
to  BT  the  difference  does  not  change.  However,  in  the  case  where  the  two 
crystals  have  different  A  constants  the  differential  frequency  change  continues 
to  increase.  This  can  be  seen  by  considering  the  difference  when  B^  =  B2, 


and  if  Fjq 


where 


AF  -  aFq  =  log(BT  +  1) 


F  A  -  F  A 
r 10  1  20  2 


is  approximately  equal  to  F2q  then, 


AF  “  aFo  =  Favg(Al  "  A2)  1o9(bt  +  1) 


F1  +F2 


avg 
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Fig.  10.  Theoretical  behavior  of  resonator  aging  using  logarithmic  rate  law 
with  different  A  constants  but  identical  B  constants  plotted  vs 
linear  time. 
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Fig.  11.  Theoretical  behavior  of  resonator  aging  using  logarithmic  rate  law 
with  different  B  constants  but  identical  A  constants  plotted  vs 
linear  time. 


Rockwell  International 


MRDC41082. 4AR 


This  expression  is  the  same  as  the  original  rate  law  with  the  A  constant 
replaced  by  the  difference  between  the  two  individual  crystal  A  constants. 

From  this  type  of  analysis  it  is  clear  that  for  a  high  degree  of  differential 
bias  stability  the  A  constants  of  the  crystals  must  be  closely  matched.  In 
fact  the  desired  result  plotted  on  a  logarithmic  time  scale  is  shown  in 
Fig.  12  and  the  frequency  change  with  increasing  time  is  zero. 

An  experimental  result  which  closely  confirmed  this  is  shown  in  the 
Appendices  of  this  report  for  the  dual  crystal  pairs  #9  &  #10,  #13  &  #14,  and 
#15  &  #16.  In  these  cases  the  individual  aging  over  30  days  was  more  than 
5  ppm,  yet  the  A  constants  of  the  two  crystals  were  sufficiently  matched  that 
the  differential  bias  stability  was  less  than  0.2  ppm  over  the  same  30  day 
period.  This  experimental  result  is  encouraging  since  it  clearly  demonstrates 
that  it  is  possible  to  fabricate  crystals  with  low  differential  bias  stability 
due  to  matched  A  constants  of  the  logramithic  rate  law. 

3.5  SAW  Sensor  Electronics  and  Inertial  Errors 

An  objective  of  this  study  was  the  determination  of  bias  stability 
errors  in  SAW  acceleration  sensors.  A  SAW  acceleration  sensor,  depicted  in 
Fig.  13,  was  designed  using  a  dual  crystal  (two  SAW  resonators  on  a  single 
crystal)  sensor.  The  dual  resonator  concept  required  two  independent 
oscillators  locked  to  the  resonator  frequencies.  The  resonator  crystal  was 
mounted  such  that  the  frequency  difference  between  the  two  oscillators  was 
proportional  to  acceleration.  The  output  from  such  a  sensor  is  inherently 
digital,  yet  quantization  errors  are  zero  since  the  counting  processes  is 
continuous.  Also  because  of  the  stability  inherent  in  crystal  controlled 
oscillators  the  dynamic  range  is  large. 

Oscillator  circuitry  for  dual  resonators  was  designed  and  constructed 
using  hybrid  chip  fabrication  techniques.  The  oscillator  circuit  design  was 
of  the  negative  resistance  type  with  one  side  of  each  resonator  returned  to 
ground.  The  negative  resistance  type  of  circuit  operated  on  the  series  resonant 
frequency  of  each  resonator  and  required  only  3  dB  of  gain.  This  was  a  direct 
result  of  using  single  port  resonators.  An  additional  advantage  of  this 
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Fig.  12.  Theoretical  behavior  of  resonator  aging  using  logarithmic  rate  law 
with  different  B  constants  but  identiacal  A  constants  plotted  vs 
logarithmic  time. 
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OSC.  #1 


13.  A  SAW  acceleration  sensor  using  dual  resonator  crystal  in 
conjunction  with  difference  frequency  measurement  system. 


Rockwell  International 


MRDC41G82. 4AR 


design  was  the  minimal  crystal  loading  which  resulted  in  high  closed  loop  Q 
and  overall  low  power  consumption.  Typically  only  10  mwatt  of  power  at  10 
volts  supply  was  required.  Also  lock-up  of  the  two  oscillators  was  not  a 
problem  even  though  the  isolation  between  them  was  only  35  dB. 

Frequency  stability  with  this  type  of  SAW  oscillator  was  quite 
good.  Shown  in  Fig.  14  is  a  typical  Allan  variance  or  fractional  deviation 
measured  using  an  HP5390  Stability  Analyzer.  Short  term  stability  was 
inversely  proportional  to  time  for  gate  times  less  than  10  msec.  For  measure¬ 
ment  times  greater  than  10  second  the  long  term  aging  and/or  temperature 
induced  frequency  fluxuations  were  noticable.  In  the  range  100  msec  to  10 
st.jnd  the  oscillators  show  a  Flicker  noise  floor  or  approximately  1x10'^. 
These  results  are  only  5  to  10  times  worse  than  the  best  bulk  wave  oscillators 
and  are  expected  to  improve. 

The  SAW  dual  crystal  oscillator  frequency  stability  determines  the 
noise  floor  or  bias  stability  of  the  SAW  sensor.  Directly  related  to  this  is 
the  sensor  instantaneous  dynamic  range.  The  full  scale  output  of  the  sensor 
is  typically  200  ppm  and  assuming  a  noise  floor  of  lxlO-10,  the  dynamic  range 
is  typically  lxlO6. 

Associated  with  this  study  was  an  evaluation  of  integrated  noise 
floors.  An  important  potential  application  of  the  SAW  sensor  is  for  inertial 
navigation  and  guidance.  Consider  the  output  of  the  sensor  being  a  train  of 
pulses  whose  count  per  second  is  proportional  tr  acceleration.  To  determine 
velocity,  V,  from  such  a  sensor  the  acceleration  count  must  be  integrated. 

This  is  equivalent  to  summing  the  total  number  of  counts  and  scaling  according 
to  the  sensor  gage  factor  and  count  duration. 

Shown  in  Fig.  15  is  the  actual  count  output  of  a  SAW  sensor  when 
measured  by  a  counter  with  a  1  second  gate  time.  The  difference  frequency  was 
60  kHz  which  has  been  subtracted  out  and  only  the  deviation  in  the  difference 
frequency  (±1  Hz)  is  shown.  Summing  the  counts  results  in  the  velocity  count 
(±10  Hz)  as  a  function  of  time  shown.  This  curve  represents  the  area  under 
the  acceleration-frequency  curve.  Performing  another  summation  as  a  function 
of  time  results  in  the  curve  for  displacement  (±4000  Hz).  The  scale  factor 
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SAW  Smaar — Duel  Oscillator 
DATE  81/11/11  TIME  11.17.00 


L4.  Allan  variance  as  measured  on  a  typical  SAW  dual  crystal  oscillator. 
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Fmmx  ■  1  Hz 

Vihx  ■  10  Hz 

Dnax  ■  4000  Hz 


Fig.  15.  Measured  count  output  of  a  SAW  dual  crystal  sensor  vs  time.  Total 
time  was  approximately  18  minutes  with  a  1  second  gate  time.  Also 
shown  is  the  computed  velocity  from  integrating  once  as  well  as  the 
displacement  as  a  result  of  integrating  the  counter  output  twice. 
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was  200  ppm/g  or  approximately  769  Hz/m/sec.2  The  integrated  velocity  error 
for  the  1000  second  time  period  shown  was  typically  less  than  0.013  m/sec  and 
the  position  error  typically  less  than  5.2  meter.  As  expected  the  integrated 
error  is  closely  related  to  the  integration  time  and  this  is  dependent  upon 
the  actual  mission  time  or  time  for  which  no  other  guidance  data  is  available. 

Sensor  integrated  errors  for  shorter  gate  and  mission  times  were 
readily  measured.  Shown  in  Fig.  16  is  a  plot  of  frequency  deviation,  velocity 
error,  and  displacement  error  when  the  gate  time  is  10  msec  for  1000  counts. 
Theoretically  this  gives  a  data  output  100  times  a  second.  The  actual  rate 
was  somewhat  slower  (0.18  sec)  because  of  dead  time  in  the  counter,  nevertheless 
the  integrated  velocity  count  was  less  than  10  Hz  and  the  displacement  count 
less  than  1000  Hz  after  a  total  integration  time  of  3  min.  This  corresponds 
to  a  velocity  error  of  less  than  0.0023  m/sec  and  a  displacement  error  of  less 
than  0.042  meters  for  a  3  min  mission  time. 

Since  the  mission  time  was  known  in  these  tests  the  cumulative  error 
per  unit  time  for  a  SAW  inertial  sensor  can  be  calculated.  In  the  first  test, 
the  mission  time  was  approximately  20  minutes,  and  the  cumulative  displacement 
error  was  5.2  meters,  giving  an  error  rate  of  approximately  15.6  meters/hour. 

In  the  second  test  the  mission  time  was  3  minutes  and  the  cumulative  position 
error  was  0.042  meters,  giving  an  error  rate  of  approximately  0.84  meters/hour. 
In  perspective  these  results  are  quite  good.  Moderately  accurate  inertial 
accelerometers  typically  have  an  error  rate  of  1  nautical  per  hour  or  1851 
meters/hour.  The  SAW  sensor  error  rates  are  considerably  below  this  value. 

Thus  the  SAW  accelerometer  noise  level  is  such  that  it  can  be  considered  as  a 
very  good  grade  of  moderately  accurate  inertial  accelerometer. 
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Fig.  16.  Measured  count  output  of  a  SAW  dual  crystal  sensor  vs  time.  Total 
time  was  approximately  3  minutes  with  a  10  msec  gate  time.  Also 
shown  is  the  computed  velocity  from  integrating  once  as  well  as  the 
displacement  as  a  result  of  integrating  the  counter  output  twice. 
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4.0  FUTURE  PLANS 

In  the  second  year  of  this  program  bias  stability,  aging,  and  SAW 
resonator  fabrication  related  studies  will  be  continued.  In  order  to  increase 
the  amount  of  data  available  construction  of  a  second  36  crystal  oven  for 
monitoring  bias  stability  has  already  begun.  In  addition  to  these  areas  of 
study  the  following  tasks  will  be  performed. 

Task  1  -  Characterize  and  mathematically  model  common  mode  rejection  of 

temperature  related  frequency  shifts  in  two-resonator  differential 
sensors.  Provided  the  structural  and  chemical  properties  of  the 
sensing  crystal  are  homogeneous,  common  mode  rejection  will  depend 
upon  mounting  stresses.  This  investigation  will  model  these  effects 
and  predict  their  magnitude  as  a  function  of  crystal  mounting  and 
proof  mass  attachment.  The  results  of  this  investigation  will  be 
coupled  with  an  investigation  of  stress  sensitivity  due  to  g-loading. 

Task  2  -  Theoretically  and  experimentally  investigate  the  sensitivity  of 
rotated  y-cut  quartz  surface  wave  resonators  to  acceleration 
effects.  The  objective  of  this  study  shall  be  to  characterize  the 
crystal  in  terms  of  second  and  third  order  piezo-elastic  constants, 
applied  proof  mass  loading,  crystal  geometry,  and  mounting  points  in 
order  to  obtain  maximum  sensitivity  to  g-loading  while  maintaining 
common  mode  rejection  of  temperature  induced  stress  effects. 

Task  3  -  To  study  long  term  bias  stability  on  prototype  SAW  accelerometers. 

This  work  will  be  an  extension  of  the  work  performed  in  the  first 
year,  however  in  this  case  the  two-resonator  crystals  will  be 
combined  with  rf  circuitry  so  as  to  provide  a  comparision  between  the 
effects  of  crystal  frequency  shifts  and  circuit  induced  frequency 
shifts.  A  goal  of  this  task  is  to  include  these  in  the  stability 
model  being  investigated. 
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-  To  study  the  acceleration  sensitivity  of  chemically  milled  surface 
wave  accelerometer  crystals.  Experimental  structures  and  fabrication 
techniques  will  be  evaluated.  Theoretical  predictions  will  be 
compared  with  experimental  results.  The  goal  of  this  study  is  to 
develop  crystal  structures  with  improved  sensitivity  to  acceleration 
while  achieving  minimal  or  zero  sensitivity  to  mounting  stresses. 

-  To  evaluate  the  overall  performance  of  prototype  SAW  accelerometers 
and  to  calibrate  the  output  in  a  standard  accelerometer  calibration 
test  system.  These  tests  will  be  performed  on  completed  sensors 
which  contain  sensing  crystals  as  well  as  hybrid  electronics. 
Calibration  of  the  sensor  will  be  performed  by  operating  the  device 
under  variable  levels  of  vibration  and  comparing  to  a  known 
accelerometer  standard.  Also  included  with  this  study  will  be  an 
evaluation  of  sensor  linearity  over  the  range  +10  g  and  a  life  test 
(fatigue)  under  vibration. 
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5.0  PUBLICATIONS  AND  PRESENTATIONS 

Two  papers  presenting  the  technical  results  of  research  performed  in 
the  course  of  this  program  were  given  during  the  last  year.  The  first  paper 
entitled  "Surface  Chemistry  Related  to  SAW  Resonator  Aging"  was  presented  at 
the  1980  Ultrasonics  Symposium  in  Boston  on  November  5-7.  This  paper 
described  the  surface  chemistry  experiments  being  performed  to  determine  which 
factors  involving  surface  chemistry  contribute  to  SAW  resonator  aging.  The 
second  paper  entitled  "Absolute  and  Differential  Aging  of  SAW  Resonator  Pairs" 
was  presented  at  the  Annual  Symposium  on  Frequency  Control.  In  this  paper  the 
results  of  our  aging  and  drive  level  studies  were  discussed  and  compared  with 
the  results  of  other  researchers  in  this  area.  For  the  purpose  of  clarity  and 
background  information,  pre-prints  of  these  papers  have  been  included  as 
appendices  in  this  report. 
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7.0  APPENDICES 

7.1  Appendix  I  -  Automatic  Aging  Measurement  Software 

This  appendix  shows  the  software  which  has  been  developed  to  auto¬ 
matically  determine  the  aging  of  quartz  resonators  using  the  specially 
designed  crystal  oven  containing  a  coaxial  switching  matrix.  A  total  of  36 
crystals  per  oven  can  be  analyzed  as  a  function  of  time.  The  data  consists  of 
a  frequency  array,  F[36,5Q],  a  string  array  for  time,  B$[50,14],  and  a  temper¬ 
ature  array,  T[50].  After  each  block  of  50  data  points  has  been  collected  the 
data  arrays  are  stored  on  magnetic  tape  (line  306).  Time  is  stored  as  a  14 
number  string  in  the  form  [MM;DD ;HH;mm;ss]  which  is  read  from  the  real  time 
system  clock  (line  112).  For  plotting  purposes  absolute  time  is  converted  to 
real  time  (line  10). 

The  data  acquisition  is  controlled  by  an  interrupt  routine  (line 
181).  How  often  the  interrupt  occurs  is  programmed  into  the  real  time  clock 
at  the  beginning  of  the  program  (line  34).  Upon  initiation  of  the  interrupt 
routine,  CYCLE  TEST  (line  95)  is  called.  This  routine  sequentially  begins  the 
testing  procedure  for  each  of  the  36  crystals.  For  each  crystal  tested  the 
frequency  of  zero  phase  (in  impedance)  is  determined  on  the  automatic  network 
analyzer.  The  temperature  is  measured  in  the  oven  and  this  is  also  stored. 

The  program  operates  from  special  function  keys  which  have  been 
previously  defined.  These  keys  allow  interactive  operation  of  the  program 
functions.  After  each  function  is  complete  the  program  asks  for  selection  of 
another  key  by  the  operator.  All  of  the  plotting  routines  for  impedance  or 
aging  data  versus  time  are  located  at  the  end  of  the  program  listing. 
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0:  "AUTOMATIC  SAW  AGING  PPOCPAM . EJS. .. 7-13-80"  : 

Is  "Software  for  36  Xtal  test  w/J. WISE's  OVENS,  trk0;file  1;  Aging  #5": 
2: 

3:  tr k  0;ldk  2 

4:  dim  F  [ 36 ,50  ]  ,  E$  [50 , 14  ]  ,T  [  50} 

5:  dim  A  $ [ 14 ]  , A [0 : 12] ,C$ [14] ,D$ [ 14] ,E$ [ 14]  ,G$ ( 32] ,H$ [36 ,2 ] , X [ 500 ] ,Y  [500] 
6:  30-A[4]  -A  (6]  -A  [  9]  *  A  [11  ]  ;28-A[2] 

7:  0-A[0]  ;31-A[1] -A [3 ) ♦ A [ 5] -A [7 ] -A [ 8] ♦A[10]  -A [12] 

8: 

9:  dsp  "SELECT  KEY  FUNCTION" ;stp  ;jmp  0 
10:  "abstime  A$":0*T 

11:  for  K=0  to  val ( A$ [ 1 ,2 ] )-l; T+A [K] *24*3600 ♦T; next  K 
12:  T+66400*val(A${4,5] )-T 

13:  T+36C0*val(A$[7,8] ) +6 C*val (A  $ [ 10 , 11 ] )+val  (A$ [ 13 , 14 ] ) -T 

14:  ret  T 

15: 

16:  "PICT  TENPErATUPE  DUPING  WAPN-UP":fxd  2;wrt  9,"F";red  9,A$ 

17:  scl  0,1000 ,20, 80;axe  0 , 20 , 100 , 1 ; axe  1000 , 8 C , 100 , 1 ; pi t  10,70,-l;lbl  A$ 

18:  for  1=1  tc  1000 

19:  ell  'MEASURE  TENPEFATUFE  ' 

20:  pit  I,X;dsp  "temperature3"  ,X 
21:  next  I; pen 

22:  dsp  "SELECT  KEY  FUNCTION" ; s tp  ;jmp  0 
23: 

24:  "PEFIODIC  DATA  ACCU ISITION" : 0-T; ent  "FFINTCUl? (YES  or  NC)",A$ 

25:  if  A$="YES" ;1*T 

26:  ell  'INTEPPUPT  SET-UP" 

27:  wrt  9, "Unit  1  Go" 

28:  ell  'DISPLAY  TINE ' 

29:  dsp  "SELECT  KEY  FUNCTION" ; stp  ;jmp  0 
30: 

31:  "DISPLAY  TINE" : 

32:  wrt  9,"Reauest  time";red  9,A$;dsp  A$;jmp  0 
33: 

34:  "INTEPPUPT  SET-UP": 

35:  wrt  9, "A"; wrt  9,"E" 

36:  wrt  9, "Unit  1=  Output  1" 

37:  er.t  "Period  (min)  ?"  ,  P;  P*6C000-*-P 
38:  fxd  0 ; wr t  9, "Unit  1  Per iod" , P, "/" 

39:  oni  9 , "INTEFPUFT" ;eir  9 

40:  ret 

41: 

42:  "NEASUFE  T  ENPEFATUPE"  :  rem  722;wrt  722  ,  "I11P6F4M3"  ;  wr  t  709,"C00E" 

43:  wait  2000 

44:  red  722,R;wrt  722 , "FlHlF2M3";wrt  709 , "03E" ; wait  100 
45:  red  722  S*lcl  722 

46-  5041.6/(ln(R)+7.15) -314 . 052-T; -3. 9526 e-7 *C ; 3 . 9747e-5-E; 1 . 2737e-8-A 

47:  C+ET  +AT*‘2  *V; S+V*X ; C-X+C 

48:  -E/2A+V(B‘2-4AC)/2A*X 

49:  ret 

50: 
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51: 

52:  "PPINT  DATA  ARRAYS" :ent  "Startinq  device  #",N 
53:  N*A;N  +  l*B;N+2-*C;  N+3-*D;N  +  4*E;  N  +  5*F 
54:  fmt  38x , "CRYSTAL  FREQUENCIES" ;wr t  701 
55:  fmt  14 x,72 ;wr t  701 

56:  fmt  4x, "TIME" , 12x, "# " , f . 0 f10x , "#" , f . 0 ,10x , M# " , f . 0 ,10 x , z 
57:  wrt  701,A,B,C 

58:  fmt  "# " , f . 0 , 10 x, "# " , f . 0 , lOx , "#" , f  .0 ,/ 

59:  wrt  701,D,E,F 

60:  fmt  86"-"-,/;wrt  701 

61:  for  1=1  tc  r0;if  I=36;dsp  "change  paper " ;beep ;s tp 
62:  fmt  cl4 ,6 fl2 . 6 ; if  I=72;dsp  "change  paper " ;beep ;s tp 
63:  if  I=r0;jmp  3 

64:  wrt  701 , E$ [ I ] ,F [A, H , F [B f I ] , F (C, 1 1 ,F [Df 1 1 , F [E , I J , F  [F, I ] 

65:  next  I 

66:  dsp  "SELECT  KEY  FUNCTION" ; stp  ;jmp  0 
67: 

66 : 

69:  "INITIALIZE  FREQ/DATA  ARRAYS": 

70:  ent  "NEW  or  OLD  data  array  to  start?",A$ 

71:  if  A$= "OLD" ;g to  "INITIALIZE  FROM  FILES" 

72:  ent  "FILE  TO  START  STORING  DATA?"  ,W 
73:  2*r0 

74:  for  1=1  to  36 

75:  fxd  0;dsp  "enter  ch#",I,"  frequency " ;ent  "",F 

76:  F/le6 -*-F  [  I ,  r0-l]  ;  next  I 

77:  dsp  "SELECT  KEY  FUNCTION" ; stp  ;jmp  0 

78: 

79:  "INITIALIZE  FFCN  FILES": 

80:  ent  "FILE  #  FOR  OLD  DATA  START?", W 
81:  rew;trk  l;fdf  W;ldf  W, F [*1 ,E$  ,T  [*] 

82:  for  1=1  to  50 

83:  if  F (1 , I ] =0 ; I+l+rO ; jrop  2 

84:  next  I 

85:  if  1  =  51 ;pr  t  "FULL  FILE";cll  'FULL  FILE  ' 

86:  ent  "PRINTOUT  STARTING  FP.EC's??",A$ 

87:  if  A$#"YES" ; jmp  6 

68:  fxd  0 ; p r  t  "STARTING  INFO:";prt  " - " 

89:  prt  "r0=  ",r0;spc 
90:  for  1=1  to  36 

91:  fxd  0;prt  "1=  ",I;fxd  6;prt  F(I ,r0-2] ♦F[I ,r0-l]  ;prt  E$[r0-l];spc 
92:  next  I 

93:  dsp  "SELECT  KEY  FUNCTION" ; beep ; s tp  ;jmp  0 
94: 

95:  "CYCLE  TEST": 

96:  for  J=1  to  36 

97:  "SET  TEST  TO  CHANNEL  l":rem  720;wrt  720, "1" 

98:  "SET  PROCESSOR  TO  C1&C2  A/R" :rem  716;wrt  716 , "Cl I4C2 I4E" 

99:  "SELECT  CRYSTAL  TO  TESTl":fmt  f.0,f.O,"E" 

100:  wrt  709 ,2 9+J-6* int ( ( J-l) /6) , 40  +  int ( ( J-l) /6) 
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101: 

102:  "INITIALIZE  FREQUENCY" : F [J , rO-1) *le6-F 

103:  ell  ' DETERMINE  ZEFO  PHASE  FREQUENCY  ' ; F/le6 *F  [ J , rO ] 

104:  "CLEAR  SCANNER" :wrt  709,"C" 

105:  next  J 

1C6:  ell  "MEASURE  TEMPERATURE" ; X-T [rO] 

107:  ell  "READ  TIME' 

108:  if  T=1 ; spc  ;spc 
109:  r  0+l-*r0  ;  wr  t  709, "C" 

110:  ret 
111: 

112:  "FEAD  TIME": 

113:  wrt  9, "Request  time" 

114:  red  9  , A$ ;  A$-*-B$  [r 0  1 

115:  ret 

116: 

117:  "FULL  FILE"  :W  +  l-*W;prt  "START  MEW  FILE";spc  ;prt  "FILE  #  =  "  ,W 

118:  for  1  =  1  to  36 ; F  II , 50 ] ♦F [1 , 1 ] ; next  I 

119:  for  1=2  to  50;for  J=1  to  3 6 ; G  +  F  [  J , I ] ; nex t  J;next  I 

120:  for  1=1  to  50 ; " " -E$ ( I ] ; Ovr  [ I ] ; next  1 

121:  2-»r0;ret 

122: 

123:  "SELECT  FESCNATCR  TO  TEST": 

124:  ent  "enter  #  to  test  1-36", N 
125:  wrt  709, "C" 

126:  fmt  1 ,2 f . 0 , "E" ;wr t  7 09 . 1  ,N+29-6* int ( (N-l) /6)  ,40+ int ( (N-l) /6 ) ; wa i t  1000 
127:  fxd  0 ; dsp  "Cev#" ,M , " :SELECT  KEY  FUNCTION" ; stp  ;jmp  0 
128: 

129:  "DETERMINE  ZEFO  PHASE  FrECUENCY": 

130:  F/le9-X 

131:  "SET  INSTRUMENTS  TO  F EMOTE" : r em  716;rem  719;rem  718;wait  1C00 
132:  "SYNTHESIZER  SET-UP":fmt  1 , , f . 0 , "  ( " ,  "4 0$ " ,  "0  0  % " ,  "500C " 

133:  "SOUPCE  CONVERTER  SET-UP" :fmt  2 , "E3M4W5S5T108V99EFA"  ,  f . 0  ,  "E" 

134:  " SOUFCE/CCNV EFTER  FREQUENCY  SET":fmt  3 , "FA" , f . 0 , "E" 

135:  "SYNTHESIZER  +  FREQ  INCREMENT" : fmt  4 , "/" , f . 0 , "A" 

136:  "SYNTHESIZER  -  FFEQ  INCREMENT" : fmt  5 , "/" , f . 0 , "B" 

137:  "SYNTHESIZER  FRECUENCY  SET":fmt  6 , "/" , f . 0 , " ( " 

138: 

139:  "SETTING  UF  SYSTEM" :dsp  "Setting  up  system" ?beeo 
140:  "SET  UP  SYNTHES IZ EF" :wr t  718.1,'f(X)' 

141:  "SET  UP  SOUFCE/CGNVERTEF" :wrt  719. 2  ,F/le6-l j ,  ait  1000 
142:  wrt  719 . 2 , F/le6+3 ; wait  1000 

143:  "ZERO  OUT  Co  REACTANCE" : dsp  "Zeroing  out  Co";beep 
144:  "F  +  lMHz*+  SYNTHESIZER  &  SOURCE" : (F+le6) /le9*X 
145:  wrt  718 . 1,  'f (X)  ';wr t  719 . 3 , ( F+le6) /le6-l ; wa i t  1000 
146:  wrt  719 . 3 , (F+le6) /le6+2 ; wa it  1000 

147:  "SET-UP  CH  1  FOR  MAG" : wr  t  716 , "C1M2S5EC1D2E" ;wait  1500 
148:  "SET-UP  CHI  ZERO  RET  LOSS":wrt  716 , "C1D4EC1D2E" ;wait  1500 
149:  "SET-UP  DISPLAY  PROCESSOR"  :wr  t  716  ,  "C2M.3S8EC2D2E"  ;wa  it  1500 
150:  "ZERO  PHASE  OFFSET" :wrt  716 , "C2D4EC2D2E" ;wait  2000 
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151:  "READ  AMP  &  PHASE":fmt  ;red  716, A, B 
152:  "CHECK  IF  PHASE  =0":if  E#0;jmp  -2 
153:  dsp  "phase  set  =" ,E;beep;wait  2000 
154:  dsp  "Resetting  start  f requency" ;beep 

155:  " RESET  START  FREQUENCY" : wr t  71 9 . 2 , F/le6-2 ; wa i t  1000;wrt  719 . 2 , F/le6+2 
156:  F/le9+X;wrt  718 .6 , 'f (X)  '; wait  1000 
157:  red  719,A;if  abs (A-F) >5e5 ; jmp  -2 
158: 

159:  "EVALUATE  PHASE  SLOPE" :dsp  "Evaluating  phase  slope" 

160:  "SET  FI— ►SYNTHESIZER" : F* 2e-6/le9*X ; wr t  718 . 5 , 'f (X)  ' ; wa it  1500 
161:  "READ  PHASE  Pl":red  716,A,B;cll  'Z-parm  phase ' ( A , P , r 1 ) 

162:  "SET  F2 ♦♦SYNTHES I Z ER" : 2X*X ; wr  t  718 . 4  , 'f (X)  '; wait  1500 
163:  "READ  FHASE  P2":red  716,A,B;cll  'Z-parm  phase' (A,E, r2) 

164:  if  r 2=rl ;gto  "RESET  START  FREQUENCY" 

165:  "CALC.  PHASE  SLOPE" :( r 2-rl) /( F* 4e-6 ) *S 

166:  fxd  2;dsp  "phase  slope  =",1000S,"  deg/kHz" ;wait  2000 

167:  if  S<0 ; gto  "DETERMINE  ZEFO  PHASE  FREQUENCY" 

168: 

169:  "ITERATE  TO  ZERO  PHASE" :dsp  "Iterating  to  zero  phase" 

170:  "SET  UP  DO  LOOP":for  1=1  to  20 

171:  "SET  NEW  FREQUENCY" :F/le9*X;wrt  718.6  ,'f  (Xf'jwait  1000 

172:  "READ  AMPLITUDE  AND  PHASE" :red  716,A,E;cll  'Z-parm  phase ' (A ,B,G) 

173:  "CHECK  IF  PHA5E=  0":if  G=0;jmp  3 
174:  "CALC.  NEW  FREQUENCY" : F-G/S*F 
175:  "END  DO  LOO P":next  I 

176:  fxd  6;dso  "F=" , F/le6 , "MHz  ", "Phase*  "  ,G , "degrees"  ;wait  1000 

177:  if  T=l;prt  "F=" , F/le6 ; pr t  "Phase*" ,G ; spc 

178:  1 c 1  7 

179:  ret 

180: 

181:  "INIEFFUPT" :beep 
182:  ell  'CYCLE  TEST' 

183:  if  r 0=51 ;cll  'RECOPD  DATA' 

184:  eir  9 
185:  iret 
186: 

187:  " f (X) "  : 

188:  1*E ;abs (X) *Y ; 0*Z 

189:  Z+Eint  (Y)  ♦  Z 

190:  10 ( Y-int ( Y) ) ♦ Y 

191:  10E*E 

192:  if  Y#0;jmp  -3 

193:  ret  Z 

194: 

195:  "PLOT  DATA  FILES" :ent  "First  file  to  plot?",P 

196:  ent  "Last  file  to  plot?",Q 

197:  ent  "Multiple  xtals  to  plot  (YES/NC) ?" ,A$ 

198:  ent  "Multiple  duals  to  plot  (YES/NO) ?" ,C$ 

199:  ent  "Multiple  diff'  to  plot  (YES/NO) ?", E$ 

200:  if  A$="YES";gto  "MULTIPLE  XTAL  PLOTS" 
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201:  if  C$="YES";gto  "MULTIPLE  DUAL  CPYSTAL  PLOTS" 

202:  if  E$="YES";gto  "MULTIPLE  DIFF  CPYSTAL  PLOTS" 

203:  "NEW  XTAL" : 

204:  ent  "crystal#(or  Temp=0)  to  plot?",X;if  X=0;cll  'PLOT  TEMPERATURE' 

205:  ent  "PLOT  SET-UP (YES  or  NO)?",A$;if  A$="NO";gto  "PLOT" 

206:  ell  'PLOT  SET-UP' 

207: 

208:  "PLOT" :dsp  "change  pen  color ?" ;beep;stp  ;if  E$="LIN" ; lire  0,E-A,R,M 
209:  if  E$="LOG" ; 1 im  0 , leg ( ( E-A ) / { 60*C ) +1) , R ,M 
210:  for  U=P  to  Q;trk  l;fdf  U;ldf  U , F [ * ]  ,B$  ,T  [  * ] 

211:  for  I=J  to  50;if  E$ ( I ] =" " ; pen;gto  216 
212:  B$  [  I  ]  -*-A$;  'abstime  A$'*Y 

213:  if  E$="LIN";plt  Y-A , ( F [X, I ] -r (X+4) ) le6/r (X+4) 

214:  if  E$*"LOG";plt  log ( ( Y-A ) / (60*C ) +1)  , ( F [X , 1 1 -r (X+4 ) ) le6/r ( X+4) 

215:  dsp  "I=",I 
216:  next  I;pen 
217:  dsp  "File=",U 
218:  next  U;lim 

219:  ent  "PLOT  LABELINC ( YES  or  NO) ?" ,A$; if  A$  =  " YES " ; 1-R ; ell  'PLOT  LABELING' 
220:  ent  "ANOTHER  XTAL  (YES/hO)  ?"  ,A$ 

221:  if  A$=" YES" ;gto  "NEW  XTAL" 

222:  dsp  "SELECT  KEY  FUNCTION" ; beep ; s tp  ;jmp  0 
223: 

224:  "PLOT  LABELINC"  : i f  E$  =  "LOG";cll  'LOG  PLOT  LABELING '; ret 
225:  csiz  3,2 

226:  if  X=0 ; pi t  0,R+(M-R)/2 ,-l;cplt  -6,0;csiz  8 ,5 , 1 ,90 ;cplt  -15/2,0 

227:  if  X=0;lbl  "Temperature  (C)”;csiz  8,5,l;gto  234 

228:  pit  (E-A)/2,0,-l;cplt  -10,-3 

229:  lbl  "1  Week  per  Division"  , 

230:  pit  0 ,R+(M-R)/2,-l;cplt  -6,0 

231:  csiz  3 ,2 ,1 ,90 ;cplt  -4,0 

232:  wrt  705,  "UC99  ,2  ,6  ,2 ,-6 ,-4  ,0,-99" 

233:  lbl  "F  (ppm)";csiz  4,2 

234:  pit  . 2 ( E-A) /2 ,M,-1 ;cpl t  0,-2;if  X=0;cplt  0,1 

235:  if  flg8;ret 

236:  ent  "plot  label?" ,G$ 

237:  lbl  C-$ 

238:  ret 
239: 

240:  "PLOT  SET-UP":trk  l;fdf  P;ldf  P, F ( *] , E$ ,T [ * ] 

241:  if  flg7 ;gto  265 
242:  1*J ;pr  t  "START  TIME®" 

243:  if  B$[l]="";prt  E$ [2] ; 2+J ;spc  ;jmp  2 
244  :  pr  t  E$ ( 1] ; spe 

245:  50-r0;trk  l;fdf  C;ldf  C»F [*J ,B$,Tl*] 

246:  pr  t  "STOP  TIME®" 

247:  if  E$(r0]  ="";r0-ln0;  jmp  0 

248:  pr t  B$ [rO] -E$;spc  ;spc  ;trk  l;fdf  P;ldf  P, F [* ] , E$ ,T [ * J 
249: 

250:  ent  "CHANGE  STAFT  TIME (YES/NO) ?" ,A$ 
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251:  if  A$="YES";ent  "ENTER  NEW  TIME (MM : DC : HH : MM : £S ) " ,C$ 

252:  if  A$="NO";B$ [J]-C$ 

253:  C$*A$; 'abstime  A$'-A 

254:  ent  "CHANGE  STOP  TIME { YES/NO) ?" ,A$ 

255:  if  A$="YES" ;ent  "ENTER  NEW  TIKE (MM: CD: HH :MM: ES ) " , D$ 

256:  if  A$="NO";E$-D$ 

257:  D$-»A$;  'abstime  A$'*B 

258:  sfg  7 

259: 

260:  ent  "maximum  +  devia  tion?"  ,M 
261:  ent  "minimum  -  deviation?" ,R 
262:  ent  "Y-tic  interval?", D 

263:  ent  "LIN  or  LOG  (T+l) ? "  ,E$ ; ent  "Time  Units (min) ?", C 

264:  dsp  "set  up  plotter " ;beep  ;stp 

265: 

266:  pclr  ;wr  t  705, "VS3" 

267:  if  flgljwrt  705, "IP1100 ,7700  ,3600 , 9700" ;cfg  l;ato  273 

268:  if  f lg2 ; wr  t  705 , " I P46 00 , 7700  ,71 00 , 9700 " ;cf a  2;gto  273 

269:  if  f lg3 ; wr  t  705 , " IP1100  ,4 700  , 3600 , 6700  "  ; cf g  3;gto  273 

270:  if  £lg4;wrt  705 ,  " IP4600 , 4 700  ,7100 , 6700 " ; cf g  4;qto  273 

271:  if  f 1 g  5 ; w  r  t  7 05 , " IP1100 , 1700 , 36 00 , 3700 " ; cf g  5;gto  273 
272:  if  f lg6 ; wr  t  7 05 , " IP46 00 , 1 700 , 7100 , 3700 " ; cf g  6 
273:  if  E$="LCC" ;gto  "LCC  SET-UP" 

274:  scl  0,E-A,R,M 

275:  wr  t  705 , "TLl , 1" ; xax  C,60*C,0,B-A 
276:  wr t  705 , "TL0 ,1 " ; xax  M,6C*C,0,E-A 
277:  if  X=0;fxd  l;csiz  5,5;jmp  2 
278:  fxd  0;csiz  2,2 

279:  yax  E-A,D,R,K;wrt  70 5 , "TLl  ,0" ; xax  R ,C*60 , 0 , E-A ; yax  C,C,F,N,1 
280:  "DETEFMINE  REF E FENCE  FFECUENCIES " : 

281:  for  1=1  to  36 ; F [I , J] ♦r ( 1+4) ; next  I 

282:  ret 

283: 

284:  "LCC  SET-UP": 

285:  scl  0 , log ( (3-A) / (60*C) +1)  ,R,M 

2b6:  wrt  705, "TLl,l" ;xax  0 , 1 ,0 , log ( (B-A) / ( 60*C ) +1) 

287:  wrt  705 , "TL0 , 1 " ; xa X  M, 1 , 0 , log ( (B-A) / ( 60*C ) +1) 

288 :  fxd  0 ;cs iz  2,2 

289:  yax  log ( ( B-A) / ( 60*C ) +1) ,C , P ,M 

290:  wrt  705 ,  "TLl ,  0  "  ;  xa  x  F.,1,0  ,log  ( (E-A)  /  ( 6C*C ) +1)  ;yax  0,D,F,M,1 
291:  gtc  "DETEFMINE.  FE FEFENCE  FFECUENCIES" 

292: 

293: 

294:  "PLOT  TEMPEFATURE" :ent  "REF  TEMP???",r99 

295:  ent  "PLCT  SET-UF  (YES/NC) ?" ,A$; if  A $= " YES" ;cll  ' 'PLCT  SET-UP' 
296:  for  U=P  to  C;trk  l;fdf  U;ldf  U,F (*] ,B$ ,T (*] 

297:  dsp  "change  pen  color ?"; beep ; s tD 
298:  for  I=J  to  50;if  E$ (I ] ="" ;pen; jmp  3 
299:  B$[IJ-A$; 'abstime  A$'*Y 
300:  pit  Y-A ,T ( I ] -r99 
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301:  next  I;pen 
302:  next  U 

303:  ent  "PLOT  LABEL ING ( YES  or  NO ) ? "  ,  A  $  ;  i  f  A$= " YES " ; gto  "PLOT  LAEELING" 

304:  gtc  204 

305: 

306  :  "RECORD  DATA": 

307:  if  W=l;rew 

308:  trk  l;fdf  Vv;rcf  W,  F  ( *]  ,  B$  ,T  [  *) 

309:  W+1*W 

310:  for  1=1  to  36 

311:  F  [  1 , 50  ]  ♦F  [1 , 1  ] 

312:  for  J=2  to  50 
313:  0*F[I,J] 

314:  next  J 
315:  next  I 

316:  T [ 50] -T[l) ;B$ [50] -E$ [1] 

317:  for  1=2  to  50 
318:  0*T[I];""*E$lI] 

319:  next  I 
320:  2-rO 
321:  ret 
322: 

323:  "1EFMINATE  PROGRAM": 

324:  wrt  9, "Abort  tiirers" 

325:  fxd  0;prt  "STOP ING  DATA" ; pr  t  "OK  FILE#  "  ,W;spc 

326:  trk  l;fdf  W;rcf  W , F  [*] ,E$  ,T  (*1 

327:  rew 

328:  wait  5000 

329:  dsp  "FEMCVE  TAPE  FROM  CALCULATOF "  ; beep  ; E  tp  ;jrrp  0 
330: 

331:  "cxd  iv”  :p3'"2+p4  “2-»p7  ;  (plp3+p2c4)  /p7+p5  ;  (p2  p3-plp4 )  /p7*p6  ;  ret 
332: 

333:  "Z-parir  phase": 

334:  tn“(pl/20) *p4 

335:  p4ccs (p2) ♦C;p4sin (p2) ♦D 

336:  ell  'cxd iv ' ( 1+C ,D,1-C,-D,A,E) 

337:  atn(E/A)+p3 

338:  ret 

339: 

340:  "LOG  PLCT  LAEELING" :cs iz  3,2 

341:  pit  log((B-A)/(60*C)+l)/2,0,-l;cplt  -11,-3 

342:  lfcl  "Logrithmic  time,  hours" 

343:  pit  0,R+(M-R)/2,-l;cplt  -6,0 

344:  csiz  3 ,2 , 1 ,90 ;cpl t  -4,0 

345:  wrt  705, "UC99 ,2 ,6 ,2 ,-6 ,-4 ,0,-99" 

346:  lbl  "F  (ppm)";csiz  4,2 

347:  pit  . 1* log ( (B-A) / ( 60*C ) +1) ,M,-l;cplt  0,-2 

348:  if  flgl  or  flg2  cr  flg3  or  flg4  cr  flg5  or  flg6  or  flg8;ret 
349:  ent  "plot  label?" ,G$ 

350:  lbl  G$ 
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351:  ret 
352: 

353:  "MULTIPLE  XTAL  PLOTS": 

354:  ent  "PLOT  SET-UP  (YES/MO) ?" ,A $ 

355:  if  A$="YES" ;cll  PLOT  SET-UP' 

356:  if  E$="LOG"  ;1  im  0 ,  log  ( (B-A )  /  (60*C ) +1)  ,R  ,M;  jirp  2 
357:  lim  0,E-A,R,M 

358:  ent  "Number  of  xtals  to  plot " , N 

359:  for  1=1  to  N;ent  "Xtal  #  ?”  ,H$ [ I ) ; next  I 

360:  for  U=P  to  C;trk  l;ldf  U , F { * J ,B $  ,T [ * ] 

361:  for  2=1  to  N 
362:  val(H$[Z])-X 

363:  for  I=J  to  50;if  E$ ( I ] =" " ; pen ; gto  368 
364:  B$[I]*A$; 'abstime  A$'*Y 

365:  if  E$="LIN" ;plt  Y-A, (F [X,I}-r (X+4) ) le6/r (X+4) 

366:  if  E$="LOG" ;  pi  t  leg ( ( Y-A ) / ( 6 C*C ) +1) , ( F [ X , I ] -r (X+4 ) ) le6/r (X+4) 

367:  dsp  "I=",I 
368:  next  I;  pen 
369:  next  Z 
370:  next  U 
371:  liir 

372:  ent  "PLOT  LABELING  (YES/NC) ?" ,A$; if  A$=  "  YES " ; 1*R ; ell  "PLOT  LAEELING ' 

373:  dsp  "SELECT  KEY  FUNCTION" ; beep ; stp  ;jmp  0 

374: 

375:  "MULTIPLE  DUAL  CRYSTAL  PLCTS": 

376:  cfg  7;sfg  l;sfg  2;sfg  3;sfg  4;sfg  5;sfg  6;sfg  8 

377:  ent  "Number  of  xtals  to  plot (12  max)",N 

378:  for  1=1  to  N;ent  "Xtal  #  ?" ,H $ [ I ] ; next  I 

379:  for  Z=1  to  N  by  2 

380:  pen#  4;cll  'PLOT  SET-UF' 

381:  if  E$  =  "LCG" ; 1 im  0 , log ( (B-A) / (60*C ) +1)  ,P  ,N ; jmp  2 
382:  lim  0  , E-A ,  F. ,  M 
383:  pen#  3 

384:  for  U=P  to  C;trk  l;ldf  U, F [*] ,B$,T [ *) 

385:  for  G=Z  to  Z+l; val (H$ [G] ) ♦X 

386:  for  I=J  to  50;if  B$ [ I ] =" " ; pe n;gto  391 

387:  B$[I]*A$;' abstime  A$'*Y 

388:  if  E$="LIN" ;plt  Y-A, (F [X , I ) -r (X+4) ) le6/r (X+4) 

389:  if  E$="LCG";plt  1 og ( ( Y-A ) / ( 6 C*C ) +1) , ( F (X , I ) -r (X+4 ) ) le6/r (X+4) 

390:  dsp  "I=”,I 
391:  next  I;pen 
392:  next  G 
393:  next  U 
394:  lim 
395:  pen#  4 

396:  ell  "FLOT  LAEELING' 

397:  fxd  0;lbl  "#" ,val (H$ [Z ] )  ,  "  6  # "  ,  val (H$ [ Z  +  l) ) 

398:  next  Z 

399:  dsp  "SELECT  KEY  FUNCTION' ;beep?stp  ;jmp  0 
400: 
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401:  "MULTIPLE  DIFF  CRYSTAL  PLCTS": 

402:  cfg  7;sfg  l;sfg  2;sfg  3;sfg  4;sfg  5;sfg  6;sfg  8 

403:  ent  "Number  of  xtals  to  plot(12  max)",N 

404:  for  1=1  to  N;ent  "Xtal  #  ?" ,H$ [ I ] ; next  I 

405:  for  Z=1  to  N  by  2 

406:  pen#  4;cll  'PLOT  SET-UP' 

,  407:  if  E$="LOG" ; 1 im  0 , log ( ( B-A) / ( 6 0*C ) +1) , R ,M ; jmp  2 

.  408:  lim  0,B-A,R,M 

409:  pen#  3 

410:  for  U=P  to  C;trk  l;ldf  U, F [ * ] ,B $ ,T [ *] 

411:  val ( H$ [ Z ] )  -»X 

412:  for  I=J  to  50;if  B$ [ I ] =" " ; pe n; qto  420 
413:  B$ [ I ] ♦A$ ; ' abstime  A$'*Y 

414:  (F [X,I ]-r (X+4) -F [X+1,I] +r (X+5) ) le6/ ( ( r (X+4) +r (X+5) )/2) ♦ F 
415:  if  E$="LCG";gto  418 
|  416:  pit  Y-A , F 

1  417:  gto  419 

418:  pit  log  ( ( Y-A ) / (60*C ) +1)  ,  F 

419:  dsp  "I=",I 

420:  next  I;pen 

421:  next  U 

422:  lim 

423:  pen#  4 

424:  ell  "PLOT  LABELING' 

425:  fxd  0 ;lbl  "#"  ,val (H$ [Z ] )  ,  "  &  # "  ,  val (H$ [Z+l] ) 

426:  next  Z 

427:  dsp  "SELECT  KEY  FU NCTICN" ; beep ; s to  ;jmp  0 
428: 

429:  "MULTIPLE  DUAL  XTAL  PCLAP  PLOT  S-PAPM": 

430:  cfg  7;sfg  l;sfg  2;sfg  3;sfg  4;sfg  5;sfg  6;sfg  8 
431:  ent  "Number  of  Xtals  to  plot {max  12)?",N 
432:  for  1=1  to  N;ent  "xtal  # ?"  ,H$ [ I ) ; next  I 
433:  ent  "File  for  reference  f reauencies?" ,U 
434:  trk  l;ldf  U , F [ * ] , E$ ,T [ *] 

435:  dsp  "Set  up  Analyzer " ;beep ; s tp 
436: 

437:  for  E=1  to  N  by  2 

438:  pen#  4;cll  "Polar  Flot  Set-up' 

439:  pen#  3 

440:  for  G=E  to  E+l; val (H$ [C] ) -X; wr t  709  , "C” 

441:  if  G=B+l;pen#  2 

442:  fmt  1 ,2 f . 0 , " E" ;wr t  7 09 . 1 ,X+2 9-6* int ( ( X-l) /6 ) ,40+int ( ( X-l) /6 ) ; wa i t  1000 
443:  F[X,l]*le6*F;F/le9»X;rem  716;rem  719;rem  718;wait  1000 
444:  fmt  1 , "/" , f . 0 , " ( " , "40$" , "CO %" , " 500C" 

445:  fmt  2 , "R3M4W5S5T108V99EFA" , f . 0 , " E" 

446:  fmt  3 , "FA", f .0 , "E" ;fmt  4 , "/• , f .0 , "A" ; f mt  5 , "/" , f . 0 , "E" 

447:  fmt  6  ,  ,  f  .0 ,  "  (  " 

448:  wrt  7 19 . 2 , F/le6-l; wait  1000;wrt  719. 2 , F/le6+3 ; wa it  1000 
449:  (F+le6)/le9-X;wrt  718 . 1 , 'f (X)  '; wr t  719 . 3 , F/le6 ;wa it  1000 
450:  wrt  719 . 3 , ( F+le6 ) /le6+2 ; wai t  1000 
j  *24330 
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451:  wrt  716 , "C1M5S2EC1D2E" ;wait  1500 
452:  wrt  716 , "C1D4 EC1D2E" ;wait  1500 
453:  wrt  716 , "C2M6S2EC2D2E" ;wait  1500 
454:  wrt  716 , "C2D4EC2D2E" ;wait  2000 
455:  fmt  ;red  716,A,T;if  T#0;jmp  -1 
456:  ell  'Zero  Phase  slope' 

457:  wrt  719 . 2 , F/le6-2 ; wa i t  1000;wrt  719. 2 ,F/le6+2 
458:  F/le9-X;wrt  718 , 6 , 'f (X)  ' ; wai t  1000 
459:  red  719,A;if  abs ( A-F ) > 5e5 ; jmp  -2 
460:  wrt  719, "W4E" 

461:  wrt  718 ,  , "41$ " , "08%" 

462:  wrt  714 , "EF" ;wr t  714 , " P3" ; wa i t  3000 

463:  for  1=1  to  250;red  714,X[I];red  714  ,Y [ 1 1 ; next  I;lcl  714 
464:  for  1  =  1  to  250  ;plt  X[I]/2  50,Y [1 1/250 
465:  next  I;pen 
466: 

467:  csiz  3 

468:  if  G=E;plt  -1 . 24 , . 9 ,-l ; f xd  0;lbl  "# " , val (H$ [ G 1  ) 

469:  if  G=B+1 ; pi t  -1. 24 , . 9 ,-l ; cpl t  0,-l;fxd  0;lbl  "#" ,val (H$ [G] ) 
470:  next  G 
471:  next  E 
472:  pen#  0 

473:  dsp  "finished  "  ;beep  ;s tp 
474: 

475: 

476:  "Polar  Plot  Set-up": 

477:  if  f lg7 ;  gto  479 

478:  dsp  "set  up  plotter "; beep ;s tp  ;sfg  7 
479:  pclr ;wr  t  705, "VS3" 

480:  if  f lgl ; wr  t  705 , " IP1100 , 7700 , 3600 , 9700 " ; cfg  ljgto  486 

481:  if  f lg2 ; wr  t  705 , " I P46 00 , 7700  ,71 00 , 9700 " ; cf g  2;ato  486 

482:  if  f lg3 ; wr  t  705 , " I P1100 , 4 700  ,3600 , 6700 " ;cf g  3;gto  486 

483:  if  f lg4 ; wr  t  70 5 , " IP46 00 , 4 700  ,71 00 , 6700 "  ;  cf g  4;qto  486 

484:  if  f lg5 ; wr  t  70 5 , " IP1100 , 1 700  ,36 00 , 3700 " ;cf g  5;gto  486 

485:  if  f lg6 ; wr  t  70 5 ,  " IP46 00 , 1700  ,7100 , 3700 " ; cf g  6 

486:  scl  -1.25, 1.25, -1,1; csiz  2 

487:  xax  0, . 1,-1 , 1 ;yax  0, . 1 ,-1,1; rad 

488:  for  1=0  to  2tt  by  ir/100;plt  cos  ( I )  , sin ( I )  ; nex t  I 

489:  pen 

490:  ret 

491: 

492:  "Zero  Phase  slope": 

493:  wrt  718 , "/" , "41$ " , "08 % " ; wr t  719 , "W4E" ; 1+Z 
494:  wr t  716 , "C2M3S8EC3M2E" 

495:  fmt  ;red  716, A, A, A 

496:  wrt  719 ,  "FC01E"  ;wait  200  ;frrt  ;red  716,A,Y[1],A 
497:  wrt  719 , "FC99E" ;wait  200;fmt  ;red  716,A,Y[2],A 
498:  if  abs(Y[l]-Y[2] )<.2;gto  501 

499:  if  Y  [  1 J  >Y[2]  jfrrt  "C30"  ,  f  3 . 0  ,  "  E"  ;  Z-l-Z  ;wr  t  716,Z;qto  495 
500:  if  Y  ( 11  <Y  [  2]  ;  f  ir.t  "C30"  ,  f  3 . 0  ,  "  E"  ;  Z  + 1*  Z ;  wr  t  7 16  ,  Z  ;gto  4  95 
501:  wrt  719 , "FC50E" ;wait  200;wrt  716 , "C2D2D4FC2D2E" 

502:  wrt  716 , "C2M6S2E" ;wrt  718 ,"/", "40$" , "00%" ; wait  1000 

503:  ret 
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7.2  Appendix  II  -  Impedance  Characteri zation 

This  appendix  shows  example  plots  of  the  impedance  of  a  typical  group 
of  36  SAW  resonators.  The  impedance  data  is  plotted  as  scattering  parameters 
(Smith  Chart  transformation).  These  are  automatically  measured  and  plotted 
using  the  measurement  software  presented  in  Appendix  I  (program  line  429). 

For  each  measurement  the  analyzer  automatically  nulls  out  the  parasitic 
capacitance  of  the  test  fixture  as  well  as  the  resonator.  This  results  in 
resonance  characteristics  nominally  symmetric  about  zero  phase  as  shown  in  the 
following  figures.  Typically  the  SAW  resonators  possess  a  series  resistance 
of  50  ohms. 
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7.3  Appendix  III  -  Absolute  Aging  Results 

This  appendix  shows  typical  data  output  for  absolute  aging  of  a  group 
of  36  SAW  resonators.  The  data  is  plotted  as  a  linear  function  of  time  for 
all  36  crystals  (program  line  353).  The  aging  data  in  this  case  has  been 
plotted  as  18  plots,  two  resonators  per  plot.  This  indicates  that  dual  reson¬ 
ator  crystals  were  tested.  The  absolute  aging  is  also  plotted  as  a  function 
of  logarithmic  time  in  the  latter  three  figures  of  this  appendix.  This  type 
of  plot  clearly  shows  the  logarithmic  straight  line  behavior  indicative  of 
chemisorption  or  oxidation. 
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Fig.  A 1 1 1 - 1  Absolute  aging  plotted  vs  linear  time  for  crystals  #1  through 

#12. 
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Fig.  AIII-2  Absolute  aging  plotted  vs  linear  time  for  crystals  #13  through 
#24. 


56 


i  Rockwell  Inte 

MRDC41082.4AR 


Z' 

■  *  25 J.  /26 

/ 

V  T  •  IF"  •'  •  v 

r 

1  Vaak  par  Oivialan 

“TT^ 

1  yna  *  #  34 

t 

• 

t 

r 

• 

« 

4 

1  . 

*  • 

t 

#  35  &  #  36 

T _ _ _ 

t  Vaah  par  Oivialan 

-t 

-• 

1  VWi  pmr  Dlvi«l«n 

-• 

Fig.  AIII-3  Absolute  aging  plotted  vs  linear  time  for  crystal  #25  through 
#36. 
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Fig.  AIII-6  Absolute  aging  plotted  vs  logarithmic  time  for  crystals  #25 
through  #36. 
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7.4  Appendix  IV  -  Differential  Aging  Results 

This  appendix  shows  differential  aging  evaluated  for  the  group  of  36 
crystals  discussed  in  the  previous  appendices.  The  data  is  first  plotted  as  a 
linear  function  of  time  in  the  first  three  figures  and  as  a  logarithmic  func¬ 
tion  of  time  in  the  latter  three  figures.  Analysis  of  differential  aging 
(program  line  401)  indicates  that  it  is  possible  to  obtain  excellent  differ¬ 
ential  bias  stability  even  though  absolute  stability  in  each  crystal  may  not 
be  possible. 
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Fig.  AIV-3  Differential  aging  plotted  vs  linear  time  for  crystals  #25 
through  #36. 
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Fig.  A I V - 4  Differential  aging  plotted  vs  logarithmic  time  for  crystals  #1 

through  #12. 


65 


Rockwell  International 


Ufpl«hat«  Mb*  haur« 


• 

• 

t 

j  ; 

1  23  t  f  24 

•%  H 

4 

Ugrt«tat*  hnr* 

4 

4 

4 

Fig.  AIV-5  Differential  aging  plotted  vs  logarithmic  time  for  crystals  #13 
through  #24. 
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jendix  V  -  Preprints  of  Pub! ications 


In  this  appendix  preprints  of  two  papers  on  absolute  and  differential 
aging  in  SAW  resonators  are  shown.  These  papers  discuss  research  carried  out 
under  contractual  support  from  DARPA. 
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SURFACE  CHEMISTRY  RELATED  TO  SAW  RESONATOR  AGING 


J.  S.  Schoenwald,  A.  B.  Marker,  W.  w.  Ho,  J.  Wise  «nd  E.  J.  Staples 


Rockwell  International  Corporation 
Thousand  Oaks,  California  91360 


Abstract 

The  electrical  characteristics  of  surface 
acoustic  wave  devices  depend  critically  on  the 
physical  and  chemical  properties  of  the  device 
surface  as  well  as  on  the  interaction  between  the 
propagating  wave  in  the  substrate  with  its  sur¬ 
rounding  gaseous  medium.  The  gas-surface  Interac¬ 
tions  were  carried  out  In  a  dual  charter  vacuum 
apparatus  with  a  SAW  resonator  Isolated  In  each 
chanoer.  One  resonator,  kept  continually  under 
vacuum,  acted  as  a  reference,  while  the  second 
resonator  was  exposed  to  a  variety  of  gases  under 
controlled  conditions.  The  nature  of  the  chemical 
Changes  occurring  on  the  surface  were  confirmed  by 
examining  the  device  surface  by  x-ray  photoelec¬ 
tron  spectroscopy.  Experiments  have  been  per¬ 
formed  to  study  the  physical  and  chemical  absorp¬ 
tion  of  reactive  and  non-reactive  gases.  These 
experiments  are  needed  to  understand  and  model  the 
aging  mechanisms  of  SAw  resonators.  Typical 
results  will  be  presented. 


1.  Introduction 

The  electrical  characteristics  of  surface 
acoustic  wave  (SAW)  devices  depend  critically  on 
the  physical  and  chemical  properties  of  the  device 
surface  as  well  as  on  the  Interaction  between  the 
propagating  wave  In  the  substrate  with  Its  sur¬ 
rounding  gaseous  medium.  These  effects,  which 
previous  researchers  have  dealt  with  only  In  terms 
of  maximizing  oscillator  and  filter  perf ormance* 
were  studied  In  this  Investigation  to  establish 
their  relationship  to  such  gas-surface  phenomena 
as  surface  oxidation,  adsorption,  and  the  Interac¬ 
tion  of  the  electroded  metallic  surfaces  and  the 
substrate  with  reactive  and  non-reactlve  gases. 

An  elementary  theory  of  kinetic  rate  pro¬ 
cesses  Is  then  derived  for  the  case  of  chemisorp¬ 
tion,  which  appears  to  be  the  predominant  effect 
observed  In  the  fabrication  and  aging  studies  dis¬ 
cussed  In  this  paper.  From  the  resulting  deriva¬ 
tion  we  obtain  Insight  concerning  the  molecular 
process  at  the  device  surface  and  how  It  brings 
•bout  a  shift  In  frequency  over  time. 

The  results  of  several  pre-aging  studies  con¬ 
ducted  under  varying  conditions  Is  presented  In 
order  to  show  the  connection  between  aging  and  the 
presence  of  reactive  gases  and  moisture  both  during 
the  fabrication  process  and  after  hermetic  sealing. 


2.  Gas-Surface  Interaction  Studies  ! 

i 

The  gas-surface  interaction  studies  were 
carried  out  In  a  dual  chamber  vacuum  apparatus 

with  a  SAW  resonator  Isolated  in  each  chamber.  i 

The  resonators  were  matched  devices  (to  within 
0.51)  operating  at  300  MHz.  One  resonator,  kept 
continually  under  vacuum,  acted  as  a  reference, 
while  the  second  resonator  was  exposed  to  a  varie¬ 
ty  of  gases  under  controlled  conditions.  The 
nature  of  the  chemical  changes  occurring  on  the 
SAW  surface  were  confirmed  by  examining  the  device  j 

surface  by  x-ray  photoeleetron  spectroscopy  < ESC A 
or  XPS).  The  SAW  resonator  devices  were  fabri-  ] 

cated  by  first  defining  the  device  pattern  by  con¬ 
ventional  photolithography  techniques  on  ST-quartz 
and  then  forming  grooves  In  the  reflector  portions 
of  the  device  via  plasma  etching.  The  transducer 
Is  masked  during  the  plasma  etching  process.  The 
aluminum  metal Izati on  pattern  of  the  reflector 
arrays  mask  the  substrate  and  etching  of  the 
quartz  substrate  takes  place  only  in  the  gaps 
between  the  metallized  stripes.  The  resonant  fre¬ 
quency  and  0  of  the  resonators  were  measured  with 
a  HP65Q5A  dual  channel  network  analyzer,  a  HP866CC 
frequency  synthesizer  and  a  HP5345A  frequency 
counter.  The  long  term  frequency  stability  of  the 
resonators  were  found  to  be  better  than  X  ppm  with 
a  Q  of  about  10.0CC.  The  measurement  accuracy  of 
the  frequency  difference  between  the  two  resona¬ 
tors  was  found  to  be  about  100  Hz. 


3.  Results  for  Non-reactlve  Gases 

The  propagation  of  surface  acoustic  waves 
•long  the  boundary  between  the  quazrtz  surface  and 
•  gas  results  In  the  generation  of  bulk  longitudi¬ 
nal  waves  In  the  gas.  This,  In  turn,  causes  an 
Increase  In  the  attenuation  along  the  direction  of 
propagation,  wlch  a  consequent  decrease  In  the  0 
of  the  resonator.  It  Is  noteworthy  that  the  shift 
In  frequency  as  a  function  of  gas  pressure  Is 
always  positive  for  the  gases  studied,  as  shown  In 
Figure  1.  As  can  be  seen,  at  a  given  gas  density, 
the  diatomic  molecules  (N?,  0?)  give  approximately 
S  times  the  shift  as  that  due  to  the  rare  gases. 
The  effect  Is  totally  reversible,  indicating  no 
significant  chemisorption  of  the  gases.  Since 
different  frequency  shifts  are  observed  for  var¬ 
ious  gases  at  a  fixed  pressure,  the  measured 
effect  cannot  be  caused  by  mechanical  distortion 
or  compression  of  the  surface  of  the  device. 

Hence,  It  can  be  concluded  that  the  perturbation 
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of  the  resonant  frequency  Is  due  primarily  to  the 
coupling  of  the  surface  eaves  to  compression  lon¬ 
gitudinal  sound  eaves  In  the  gases.  Although 
there  Is  no  existing  theory  for  quantitatively 
describing  this  phenomena.  It  Is  expected  that  the 
magnitude  of  the  shift  Is  related  to  the  compress¬ 
ibility  and  velocity  of  the  gas  present.  It  Is 
therefore  reasonable  to  expect  that  diatomic  gases 
•111  have  a  larger  effect  than  the  rare  gases. 
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Figure  1  -  Frequency  Shift  of  St-Quartz  Resonators 
as  a  Function  of  Gas  Loading  at  300  l*iz. 


These  results  Indicate  that  great  care  must  be 
exercised  in  vacuum  or  noble  gas  encapsulation  of 
SAW  devices  since  air  leakage  can  result  In  fre¬ 
quency  shifts  as  large  as  30  ppm. 

The  attenuation  of  the  surface  acoustic  waves 
caused  by  the  emission  of  compressions!  waves  In 
the  gases  was  also  determined  by  measuring  the 
Changes  In  the  quality  factor  Q  of  the  SAW  resona¬ 
tors  as  a  function  of  gas  density.  Typical  results 
Of  the  attentuatlon  measurements  for  nitrogen  gas 
are  shown  in  Figure  2  as  a  function  of  gas  pres¬ 
sure.  The  data  Is  linear  in  pressure  and  has  a 
value  for  the  slope  of  2. OS  >  10*®  db/cm-MHz- 
■nHg.  Similar  data  were  obtained  for  02  which 
showed  a  slope  of  2.14  ■  10*®  db/cm-MHz-nriHg.  The 
measurements  showed  good  agreement  for  repeated 
cycling  of  the  gases.  Since  oxidation  of  fresh 
olumlnum  surfaces  Is  expected  to  occur  in  the 
presence  of  oxygen,  the  absence  of  any  non-reversl- 
ble  effect  demonstrates  that  the  surface  Is  In  the 
form  of  A120j  and  no  further  surface  reaction  Is 
taking  place. 

There  are  several  existing  theories  describing 
the  attenuation  loss  In  propagating  surface  waves 
along  the  boundary  of  a  solid  substrate  and  a 
gas.  These  are  the  approximation  approach  devel¬ 
oped  by  Artz  and  co-workers,*  the  continuum  mechan¬ 
ical  approach  originally  proposed  by  Campbell  and 
Jones*  and.the  perturbation  formulation  by  Tlersten 
and  Slnhe.®  The  first  two  approaches  have  been 
used  by  Slobodnlk*  for  comparison  with  experimental 
data  obtained  with  various  monatomic  gas  loading  on 
T-cut  and  Z-propagatlng  LINbOj  surfaces.  Me  found 
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Figure  2  -  Surface  Wave  Attenuation  as  a  Function 
of  W2  Gas  Loading  for  St-Quartz  SAw 
Resonator  at  300  MHZ . 


that  the  theoretical  calcualtlons  bracketed  the 
experimental  results  with  the  approximate  theory 
giving  consistently  higher  values  for  the  attenua¬ 
tion.  When  these  same  calculations  are  scaled  to 
the  case  for  02  and  W2  loading  on  ST-quarts  at 
300  MHz,  the  correct  pressure  dependence  is  ob¬ 
tained  but  the  absolute  magnitudes  are,  respective¬ 
ly,  20t  and  SOt  lower  than  the  experimental  val»es 
observed  In  the  present  Investigation. 

The  more  recent  theory  of  Tlersten  and  Sinna,5 
Include  the  previously  neglected  effect  of  gas  vis¬ 
cosity  which  significantly  increases  the  attenua¬ 
tion.  The  result  of  their  calculation  for  air 
loading  on  ST-quartz  at  300  Wz  Is,  however,  still 
about  30t  lower  than  the  present  results. 


4.  Results  for  Reactive  Gases 

Experiments  were  carried  out  to  study  the  ad¬ 
sorption  and  chemisorption  of  gases  on  the  Saw 
resonator  surface.  Since  extensive  previous  work 
has  been  carried  out  In  this  laboratory  on  the  gas- 
surface  chemistry  of  HC1  gas  and  A120i  surfaces, 
this  system  was  chosen  for  the  Initial  experi¬ 
ment.  The  test  resonator  surface  was  exposed  to 
known  concentrations  of  MCI  with  N2  as  the  carrier 
gas  In  the  measurement  system.  The  changes  In  the 
electrical  properties  of  the  resonator  were  mea¬ 
sured  In  situ  as  a  function  of  exposure  time.  The 
resonator  surface  was  chemically  characterized  by 
XR$  before  and  after  the  exposure. 

Upon  exposure  of  the  test  SAW  device  to  MCI 
fas,  the  resonant  frequency  shifts  rapidly  towards 
a  lower  frequency.  This  trend  then  continues  as  a 
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function  of  time  at  •  constantly  decreasing  rate. 
When  the  system  Is  evacuated,  the  frequency  Immedi¬ 
ately  returns  to  a  steady  state  value  which  Is  sig¬ 
nificantly  lower  than  the  Initial  vacuum  resonant 
frequency  for  the  device  before  exposure.  Upon  a 
second  exposure  of  the  device,  the  resonant  fre¬ 
quency  again  decreases  and  shows  the  general  behav¬ 
ior  as  that  during  the  first  exposure.  However, 
upon  evacuation,  the  vacuum  resonant  frequency  Is 
again  shifted  to  a  lower  value  from  that  before  the 
exposure. 

The  data  from  a  Typical  run  (It  HC1  In  No. 

25  cm  total  pressure)  are  shown  In  Figure  3,  in 
which  the  background  shift  to  higher  frequency 
caused  by  the  Nz  gas  loading  Is  corrected  for  via 
the  data  shown  previously  In  Figure  1.  As  can  be 
seen,  the  irreversible  shift  In  the  vacuum  resonant 
frequency  also  decreases  with  each  successive 
cycling  of  the  gas.  Hhen  the  concentration  of  HC1 
gas  Is  varied,  it  was  found  that  the  measured 
effect  scales  linearly  with  HC1  concentration. 
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Figure  3  •  Resonant  Vacuum  Frequency  Shift  as  a 
Function  of  Exposure  Time  for  It  MCI 
Concentration  In  Nitrogen  Sas  at  25  cm 
HB  n2  at  Room  Temperature  (296k).  Also 
Shown  Is  HC1  Adsorption  Rate  on  AUO3 
Surfaces  as  Determined  by  Net  Chemistry 
Analysis  (note  that  of  Is  negative). 


XPS  analysis  of  the  SAW  resonator  before  and 
after  MCI  exposure  showed  a  chlorine  uptake  and  a 
loss  of  oxygen  from  the  AUO3.  The  ratio  of  the 
chlorine  uptake  to  oxygen  Toss  was  2  to  1,  Indi¬ 
cating  that  the  AlyOv  was  being  chemically  con¬ 
verted  to  AlClj.  Further  confirmation  that  the 
chemisorption  was  occurring  on  the  Al^Oj  rather 
than  on  the  quartz  surface  was  provided  by  addi¬ 
tional  XPS  studies  made  on  a  blank  quartz  substrate 
before  end  efter  exposure  to  MCI  gas  under  Identi¬ 
cal  conditions.  The  results  showed  no  significant 
Cl*  present. 


The  electrical  neasureeients  together  with  the 
XPS  findings  demonstrate  that  both  reversible  ad¬ 
sorption  and  Irreversible  chemisorption  of  HC1  was 
occurring  on  the  AUO3  surface.  The  observed  time 
dependence  of  the  data  given  In  Figure  3  Indicate 
that  when  fresh  AI-O3  surfaces  are  exposed  to  HC1, 
an  Initial  gas-surface  reaction  occurs  which  forms 
a  chloride  coating.  This  process  Is  then  followed 
by  adsorption  of  gaseous  HC1  onto  bound  chloride 
ions  until  equilibrium  is  reached  between  the  vapor 
pressure  of  the  adsorbed  HC1  coating  and  the  gas 
phase  HC1  concentration.  At  this  point,  the  sur. 
face  HC1  concentration  saturates  and  no  further 
adsorption  occurs.  Upon  evacuation  to  vacuum, 
desorption  of  a  surface  takes  place  and  the  shift 
In  the  vacuum  frequency  of  the  resonator  Is  a  mea¬ 
sure  of  the  net  chemisorption  on  the  surface. 

In  order  to  quantitatively  relate  the  electri¬ 
cal  measurements  to  the  surface  chemistry,  the  data 
are  compared  to  the  results  from  a  previous  study 
undertaken  In  this  laboratory  to  determine  the  sur¬ 
face  absorption  efficiency  (the  ratio  of  reactive 
collisions  to  the  total  number  of  gas-surface  col¬ 
lisions)  of  HC1  gas  onto  alumina  particles.  In 
this  prior  work,  total  adsorbed  Cl*  on  AI-O3  parti¬ 
cles  were  measured  as  a  function  of  both  RCT  con¬ 
centration  and  exposure  time  with  a  wet  chemical 
analysis  technique,  the  data  were  found  to  be  In 
agreement  with  the  expression 


d[HCl]s/dt  •  yVA[HCl]0/4  (1) 


where  V  Is  the  mean  kinetic  velocity  of  the  HC1  gas 
molecule,  A  Is  the  A1Z03  surface  area,  [HC1].  Is 
the  vapor  phase  concentration  an C  [h:?J  is  the 
total  adsorbed  HC1  on  the  surface.  The  observed 
HC1  adsorption  as  a  function  of  time  from  this  pre¬ 
vious  work  Is  scaled  to  the  experimental  conditions 
of  the  present  study  and  the  excellent  agreement  in 
the  general  shape  of  the  curves  seen  In  Figure  3 
demonstrates  that  the  present  experimental  method 
can  be  used  for  quantitative  real-time  measurement 
of  the  adsorption  processes. 

The  sensitivity  of  the  measurement  method  can 
be  evaluated  by  a  direct  comparison  of  the  data 
given  In  Figure  3.  The  result  shows  that  the  fre¬ 
quency  shift  Is  related  to  the  absolute  BC1  adsorp¬ 
tion  via  the  calibration  factor,  3  *  10*'  mole- 
cules/cnr  ■  1  kHz,  corresponding  to  a  sensitivity 
for  HC1  adsorption  at  the  level  of  3  to  5  mono¬ 
layers  on  the  device. 

Similar  experiments  were  conducted  to  deter- 
mlne  the  frequency  shift  caused  by  exposure  of  the 
SAN  resonator  to  water  vapor  under  atmospheric  con¬ 
ditions.  Measurements  were  made  with  a  SAW  resona¬ 
tor  having  a  clean  surface  as  well  as  one  which  has 
been  previously  exposed  to  HC1  gas.  It  was  found 
that  the  lowering  of  the  SAN  frequency  was  totally 
reversible  for  the  former  case,  whereas  a  large 
Irreversible  vacuum  frequency  off-set  was  observed 
for  the  latter  case.  The  magnitude  of  the  shifts 
were  considerably  larger  than  that  found  for  HC1 
end  the  time  scale  for  the  effect  was  approximately 
a  factor  of  bro  shorter.  The  results  from  a  typl- 
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C«1  run  at  75»  relative  humidity  art  shown  In 
Figure  4.  The  total  observed  shift  was  about 
40  kHz  and  the  observed  vacuum  offset  before  and 
after  exposure  was  14  kHz.  This  sensitivity  to 
water  vapor  Is  most  likely  caused  by  H20  deposition 
over  the  entire  device  surface,  rather  than  react¬ 
ing  specifically  with  the  aluminum  as  Is  the  case 
for  NCI.  The  experimental  results  obtained  for  the 
clean  surface  Indicate  that  no  chemisorption  Is 
taking  place  for  water  vapor  on  the  quartz  or  the 
oxide-coated  aluminum  surfaces.  The  observed  chem¬ 
isorption  for  the  device  previously  exposed  to  HC1 
vapor  can  be  attributed  to  the  highly  deliquescent 
nature  of  AlClv  which  readily  forms  AlClj.Bh^O  In 
the  presence  or  water  vapor. 


Figure  4  -  Resonant  FrequencyShlft  as  a  Function 
of  Exposure  Time  for  75S  Relative 
Humidity  In  Air  at  Atmospheric  Pressure 
and  Room  Temperature  (296K). 


5.  Rate  Theory  of  Kinetic  Processes 

Parker  has  pointed  out  that  frequency  aging 
effects  due  to  chemisorption  and  oxidation. behave 
according  to  the  logarithmic  rate  relation7 


af/f  -  A*log(Bt  ♦  1)  (2) 


bq( t)  •  log(l  ♦  t/t0) 


q  Is  the  number  of  particles  chemisorbed  per  unit 
area  at  time  t  after  the  beginning  of  the  experi- 
aent.  t.  Is  a  time  constant  depending  on  the  phys¬ 
ical  conditions  and  b  Is  a  constant  appropriate  to 
the  quantity  being  Matured. 

Suppose  that  at  the  time  t  •  0  •*  have  a  sun- 
face  with  adsorption  sites  in  contact  with  a  gas. 
The  nature  of  these  sites,  and  their  node  genera¬ 
tion,  Is  not  specified  by  the  theory,  nor  is  the 
amount  of  gas  already  adsorbed  on  the  surface  at 
time  t  •  0  of  Importance.  The  sites  may  be  due  to 
the  properties  of  that  part  of  the  surface  which 
has  not  yet  been  covered  by  adsorbed  molecules, 
they  may  be  produced  by  the  action  of  the  already 
adsorbed  molecules  on  the  surface,  and  they  might 
even  be  offered  by  adsorbed  molecules  themselves. 

At  this  time  let  there  be  s0  sites  per  unit  area, 
and  let  there  be  $  sites  per  unit  area  at  an  arbi¬ 
trary  time  t.  Let  there  be  N  Impacts  by  gas  mole¬ 
cules  with  the  surface,  expressed  per  unit  area  and 
per  unit  time.  Let  a  be  the  effective  contact  area 
between  a  molecule  and  the  surface  upon  collision, 
to  that  a*s  Is  the  fraction  of  Impacts  that  could 
be  effective  in  trapping  a  molecule.  One  may  ex¬ 
pect  a  to  be  roughly  of  the  order  of  a  molecular 
cross  section.  The  probability  that  a  molecule 
sticks  upon  Impact  is  proportional  to  a*s,  i.e.. 


p  •  Cas  (4) 


where  C  Is  a  constant.  The  number  of  molecules 
Sticking  per  unit  area  per  unit  time  Is 


dq/dt  •  CasN  .  k(5) 


Let  b  be  the  effective  area  over  which  sites  become 
Invalidated  by  the  adsorption  of  a  single  mole¬ 
cule.  b  Is  not  necessarily  of  the  same  order  as  a, 
since  It  is  the  area  covered  by  an  adsorbed  mole¬ 
cule  (which  is  of  the  order  of  a)  reduced  by  the 
area  of  a  new  site  generated  on  top  of  the  adsorbed 
molecule.  This  Is  reasonable  to  expect  since  It  Is 
likely  that  a  surface  of  adsorbed  molecules  are 
likely  to  offer  new  sites  to  further  adsorption. 

The  number  of  sites  Invalidated  per  unit  area  per 
unit  time  Is  therefore 


where  A  and  B  are  fit  parameters.  Equation  (2)  Is 
a  specific  statement  of  the  more  generalized 
Clovlch  equation,  which  describes  catalytic  oxida¬ 
tion  and  chemisorption  of  various  gases.on  a  var¬ 
iety  of  aetal  and  Mtal  oxide  surfaces.8 

We  will  now  derive  the  general  form  of  the 
Clovlch  equation  in  order  to  appreciate  the  physi¬ 
cal  basis  of  the  parameters  A  anbd  B  In  (2).  We 
begin  by  first  presenting  the  more  general  fon  of 


-ds/dt  •  CNas*bs  .  (6) 

Integration  yields 

s0/s  ■  1  ♦  C*iabset  .  (7) 

Solving  for  s.  Inserting  In  (S)  and  solving  for  q 
yields 
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q  •  (l/b)1n(I  ♦  (t/te)).  1/t.  -  Ceabs. . 

(8) 

This  Is  the  logarithmic  rate  law  known  as  the 
Elovlch  aquation.  Remembering  that  b  Is  the  area 
Over  which  no  further  adsorption  can  take  place 
because  there  Is  no  catalyst,  we  may  aultiply  both 
sides  of  (8)  with  a  constant  of  our  choice  to 
express  the  dependence  of  some  quantity  on  the  rate 
of  adsorbed  molecules.  Thus,  If 


.  m  fractional  frequency  shift 
*  *  adsoroed  molecules/uni t  area 


then  A-q  •  af/f,  A'/b  ■  A  and  (l/t0)  •  8  as  defined 
by  Parker  In  Equation  (1). 


6.  Results  of  Aging  Studies 

Several  lots  of  375  MHz  single  pole  SAW  reso¬ 
nators  were  fabricated,  using  a  process  previously 
described.  The  hermetically  sealed  devices  were 
maintained  In  an  oven  at  the  turning  point  tempera¬ 
ture  corresponding  to  the  cut  of  quartz  and  metal¬ 
lization  thickness  -  typically  65-75  degrees  centi¬ 
grade.  Up  to  six  crystals  could  be  probed  via  a 
coaxial  rf  switch.  Switching  between  the  crystals 
Is  accomplished  by  an  HP3495  Scanner,  which  also 
compares  the  oven  thermocouple  to  Its  Own  Internal 
standard  to  measure  oven  temperature.  The  resona¬ 
tors  are  probed  by  Instructing  the  scanner  to  acti¬ 
vate  the  rf  switch  to  select  the  desired  crystal. 

7  •  entire  system  1*  controlled  by  an  HP9825A  sys¬ 
tem  controller  via  an  IEEE-488  data  bus.  By  compu¬ 
ter  algorithm,  the  system  controller  find  the 
series  resonant  frequency  of  each  crystal,  measures 
the  oven  temperature  and  logs  the  time  and  data, 
using  the  HP9&G35A  real  time  clock.  The  system  can 
be  programed  to  'wake  up*  on  an  Interrupt  comnand 
from  the  clock  and  complete  the  cycle  of  measure- 
arnnts  on  all  the  crystals.  An  HP8505A  network  ana¬ 
lyzer  measures  crystal  return  loss  ($},)  to  obtain 
electrical  characteristics  and  resonant  frequen¬ 
cy.  Data  Is  stored  on  magnetic  tape  and  maybe  re¬ 
trieved  at  any  time  for  plotting,  analysis  or 
further  data  reduction,  figure  5  Is  a  schematic  Of 
the  data  collection  system,  with  this  arrangement 
the  crystals  remain  undisturbed  by  temperature  var¬ 
iations  or  vagaries  of  Mnual  data  collection. 

Oata  collection  can  proceed  In  the  absence  of  any 
personnel . 

The  frequency  aging  results  of  several  repre¬ 
sentative  sets  of  SAW  resonators  will  be  discussed 
In  the  remainder  of  this  section.  In  the  first  set 
of  crystals,  a  great  deal  of  variation  in  the  fab¬ 
rication  process  took  place,  even  though  the  gen¬ 
eral  format  of  the  process  was  followed,  f.c.,  de¬ 
vices  came  from  different  wafers,  were  plasma 
etched  Individually  for  different  periods  of  time, 
were  subjected  to  repeated  wet  add  clean-up  pro¬ 
cedures,  and  were  Irregularly  treated  with  regard 
to  pre-  and  post-sealing  vacuum  bakeout.  The  frac¬ 
tional  shifts  found  to  occur  In  these  devices  Is 
shown  In  figure  6.  The  oven  temperature  was  nom¬ 
inally  75  degrees  C,  with  a  characteristic 


fluctuation  of  t0.2  degrees.  The  Interesting 
feature  of  figure  6  Is  not  the  differences,  but  the 
similarities.  All  curves  have  the  same  character¬ 
istic  shape,  except  for  some  competing  mode  of 
aging  In  two  of  the  crystals  which  damped  out 
within  the  first  24  hours.  It  was  found  that  the 
crystal  with  the  worst  aging  also  had  the  lowest  C. 
and  comparison  with  the  crystal's  performance  prior 
to  vacuum  sealing  indicates  that  the  hermetic 
Integrity  of  the  package  was  not  maintained. 

Blanks  In  the  data  collection  occur  because  the 
system  controller  was  In  use  for  other  purposes, 
but  the  trend  in  frequency  aging  Is  clearly  evident 
for  the  one  month  period  during  which  data  collec¬ 
tion  occurred. 


figure  6  -  fractional  Shift  In  frequency  of  Six  SAW 
Resonators  Aging  at  75*C  for  four  weeks. 

Mid-way  through  the  data  collection  period  the 
oven  was  cycled  off  and  on  between  continuous 
Intervals  of  measurement.  A  nonrecoverable  up¬ 
shift  occurs  In  all  crystal  frequencies  due  to  the 
temperature  retrace,  quite  apart  from  the  smooth 
aonotonic  drift  the  crystals  with  time  at  constant 
temperature.  It  Is  also  noteworthy  that  the  crys¬ 
tals  with  higher  aging  rates  also  have  higher  non¬ 
recoverable  up-shifts  due  to  temperature  retrace. 
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At  mentioned  earlier,  there  It  ttrong  evidence  thet 
this  shift  It  due  to  the  pretence  of  water  vtpor 
end  oxygen  or  other  reactive  gates  In  the  package 
due  to  lack  of  hermltlclty  or  process  control. 
Out-gastlng  of  aol sture  diffused  In  the  quartz  lat¬ 
tice  may  be  a  significant  factor.  The  correspon¬ 
dence  between  high  device  0.  low  aging,  and  very 
sanal  1  shifts  In  frequency  due  to  temperature  re¬ 
trace  seem  to  be  coupled  to  the  rate  at  which  reac¬ 
tive  chemical  processes  occur  at  the  crystl  sur¬ 
face,  and  these  Inturn  are  related  to  the  partial 
pressures  of  the  gasses  In  the  crystal  enclosure. 
Further  evidence  In  support  of  this  phenomenon  will 
be  presented  later  In  this  paper. 

The  dependence  of  aging  on  environmental  mois¬ 
ture  and  other  gasses  that  might  contribute  by  re¬ 
action  or  oxidation  at  the  crystal  surface  was 
studied  by  preparing  a  set  of  samples  In  the  usual 
manner,  with  the  following  exception.  A  24  hr  high 
vacuum  bakeout  at  ISO  degrees  C  preceded  resistance 
weld  sealing  under  high  vacuum.  The  glove  box  en¬ 
closing  the  champer  was  purged  of  moisture  to  the 
extent  that  the  dew  point  was  -40  degrees  C.  In  an 
environment  relatively  free  Of  moisture,  nonre- 
verslble  aging  effects,  such  as  those  discussed 
earlier  In  which  resonators  were  exposed  to  mois¬ 
ture  after  prior  exposure  to  NCI,  ft  Is  expected 
that  aging  rates  should  be  considerably  lower. 
Figure  7a  shows  the  results  for  a  set  of  six  crys¬ 
tals  sealed  under  high  vacuum  after  a  24  hour  bake¬ 
out  In  high  vacuum.  The  long  term  fluctuations  are 
caused  by  24  hour  cyclic  fluctuations  in  the  oven 
temperature,  as  can  be  seen  In  Figure  7b  (later 
remedied),  and  the  short  term  fluctuations  are  due 
to  the  relative  shift  of  the  turning  point  for  this 
set  of  devices  (50  degrees  C)  from  the  oven  tempea- 
ture  (73  degrees).  It  can  readily  be  seen,  how¬ 
ever,  that  w ithln  two  days,  crystal  aging  subsides 
below  0.5  ppm/day  after  temperature  fluctuation 
effects  are  normalized  out.  This  Is  an  order  of 
siagnltude  Improvement  over  aging  rates  for  previous 
devices. 

Further  confirmation  of  device  sensitivity  to 
ambient  atmospheric  conditions  was  revealed  during 
the  frequency  fine  tuning  procedure.  The  resonant 
frequency  of  a  device  Is  tracked  with  an  HP8660C 


tm  (wni 


- '■ - 1 - 1 - i _ a  a _ l 


irt 


Figure  7  -  Low  Aging  Rate  Saw  Resonators. 


synthesized  signal  generator,  an  HP840SA  and  a  dual 
directional  coupler.  A  thermocouple  beneath  the 
TO- 5  header  monitored  temperature.  After  fine 
tuning  the  device  frequency  In  the  plasma  etcher, 
the  chamber  was  repressurized  via  an  air  Inlet 
valve.  Immediate  re-evacuation  of  the  chamber, 
followed  by  a  suitable  time  delay  to  relax  tempera¬ 
ture  fluctuations  during  pressure  changes,  resulted 
In  a  nonreverslble  shift  up  In  frequency  by  30-43 
kHz  for  resonators  operating  at  375  Mhz.  This 
result  Is  consistent  with  that  obtained  for  300  Miz 
resonators  exposed  to  moisture  after  exposure  to 
reactive  gases.  No  epoxy  or  other  mounting  rwte- 
rial  was  used;  Instead,  the  devices  were  strapped 
to  the  header  using  gold  wire.  Therefore,  contam¬ 
ination  from  mounting  material  was  not  a  factor. 

The  probable  cause  of  this  effect  Is  the  removal  of 
at  least  part  of  the  oxidized  or  chemisorbed  layer 
from  the  aluminized  portions  of  the  device  during 
plasma  etching  In  CF.,  which  thins  the  metal  film 
(raising  velocity  and  therefore  frequency), 
followed  by  reoxidation  in  the  presence  of  atmo¬ 
spheric  humidity  (further  raising  frequency). 

A  frequency  shift  of  this  magnitude  (100  ppm) 
Is  catastrophic  In  a  device  where  the  Q  Is  typi¬ 
cally  greater  than  10,000.  Clearly,  the  presence 
of  moisture  and  reactive  gases  at  the  crystal  sur¬ 
face  during  the  course  of  fabrication  and  sealing 
plays  a  demonstrated  role  In  the  absolute  frequency 
and  long  term  stability  of  SAW  resonators.  Envi¬ 
ronmental  conditions,  certainly  in  the  final  stages 
of  fabricating  a  tuned  sealed  device,  should  be 
carefully  controlled. 

Batch  handling  of  crystals  through  the  entire 
fabrication  cycle  except  the  final  tuning  provides 
additional  evidence  that  uniformity  of  the  envlron- 
mental  history  for  all  devices  results  in  uniform¬ 
ity  of  long  term  stability.  Figure  8  shows  aging 
results  for  a  group  of  crystals  we  have  recently 
begun  to  Investigate.  The  entire  group  was  jfro- 
cesses  as  a  single  wafer  section  and  segmented  Into 
chips  only  after  the  Initial  plasma  etch  and  just 
prior  to  mounting.  No  extraordinary  procedures 
were  undertaken,  such  as  long  term  high  vacuum 
backout.  Only  In  the  final  tuning  and  sealing 
steps  were  crystals  treated  Individually,  but  the 
procedure  for  each  was  relatively  uniform.  The 
resultant  aging  rates  are  very  tightly  clustered. 
Low  Initial  aging  rates  was  not  the  objective  in 
this  study.  Rather,  we  were  concerned  with  the 
variation  that  might  be  expected  In  a  very  uniform¬ 
ly  treated  set  of  samples,  and  we  observe  that.  In 
fact,  the  aging  Is  quite  uniform. 


Figure  0  -  Preaging  of  Six  Crystals  Subjected  to 
Very  Uniform  Process  Conditions. 
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A  comparison  with  the  aging  Equation  (1)  for 
A  •  3  ppm  and  1.5<6<3.5  plotted  In  figure  9  shows 
how  closely  the  aging  data  In  figure  6  follows  a 
logarithmic  rate  dependence  with  t0  ■  1/8  «  20-30 
ailnutes  being  the  time  constant  character) stlc  of 
the  chemlsorptlve  process  occurring  at  the  elec¬ 
trode  surface. 


figure  9  -  A  ln(Bt  ♦  1)  Plotted  for  Values  of  A,  8 
that  Model  Aging  Results  Shown  In 
figure  8. 


7.  Conclusions 

He  have  presented  results  of  the  effects  of 
non-reactive  and  reactive  gases  on  the  frequency 
and  attenuation  of  SAW  resonators,  for  non¬ 
reactive  gases  It  has  been  determined  that  the 
effect  of  the  gas  Is  totally  reversible  upon  return 
of  the  SAW  resonator  to  a  condition  of  high 
vacuum.  If  complete  oxidation  of  the  aluminum 
electrodes  has  preceded  measurement,  the  effect  of 
oxygen  is  also  completely  reversible.  As  an  ex¬ 
ample  of  reactive  gas  effects,  exposure  of  a  SAW 
resonator  to  HC1  gas  which  had  previously  been  ex¬ 
posed  to  oxygen  to  allow  electrode  oxidation  re¬ 
sulted  In  a  nonreverslble  frequency  shift  down. 
Subsequent  exposure  to  moisture  In  air  results  In  a 
nonreversible  shift  up.  The  Initial  exposure  to 
MCI  converts  AKOj  to  AlClj.  Moisture  then  com¬ 
bines  to  form  AICT3-6H2O.  X-ray  photoelectron 
spectroscopy  (XPS)  was  used  to  confirm  the  presence 
of  the  various  compounds  mentioned.  Comparison 
with  blank  substrates  of  quartz  Indicate  that  no 
chemical  process  occurs  In  which  chlorine  Is  taken 
up,  from  which  we  conclude  that  the  aluminum  1$  the 
catalytic  surface  promoting  the  chemi sorptive  reac¬ 
tion.  Wet  chemistry  analysis  Indicates  that  the 
observed  frequency  shifts  result  from  the  chemi- 
sorptive  formation  of  3  to  5  monolayers  on  the 
device. 

We  have  presented  an  elementary  theory  of 
kinetic  rate  processes  as  applied  to  chemisorp¬ 
tion.  This  results  In  a  molecular  level  model  that 
predicts  the  logarithmic  aging  rate  behavior  of  the 
SAW  resonators. 

We  have  presented  the  results  of  SAW  resonator 
pre-aging  studies  under  various  conditions  and  have 
shown  the  connection  between  environmental  control 
of  the  fabrication  process  and  device  stability 
after  fabrication.  The  effects  of  moisture,  cata¬ 
lytic  oxidation  and  reactive  gases  can  be  measured 
and  quantified.  It  Is  shewn  that  steps  can  be 
taken  to  reduce  the  effect  of  these  aging  mechan¬ 
isms  on  the  long  term  stability  of  SAW  resonators. 
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ABSOLUTE  AND  DIFFERENTIAL  AGING  OF  SAW  RESONATOR  PAIRS* 
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This  paper  examines  the  absolute  and  differ¬ 
encial  aging  of  SAW  resonators  fabricated  as 
pairs.  Low  absolute  aging  rates  of  individual  SAW 
resonators  have  been  observed  after  hermetic  seal¬ 
ing  and  temperature  stab  llzacion  inside  70°C  pre¬ 
aging  ovens.  Typical  drift  rates  fall  below 
0.5  ppa/day  within  48  hours  and  0.1-0. 5  ppm/wk 
afcer  3  weeks.  The  differential  drift  rate 
between  pairs  was  less  than  0.25  ppm/day  afcer  48 
hours  and  typically  0. 1-0.2  ppm.'wk  after  3  weeks. 

The  effect  of  drive  level  was  studied  for 
devices  fabricated  using  evaporated  chromium 
( 50i)/alumlnum  (IOOOA),  rf  sputtered  aluminum 
(1000S)  and  evaporated  aluminum  +  4 . 5 "X  copper 
alloy.  Significant  performance  degradation  occurs 
at  75°C  for  power  levels  above  +10  dBm  delivered 
to  the  crystal  at  resonance.  Comparative  studies 
of  frequency  dependence  on  drive  level  show  no 
significant  differences  among  the  three  metal¬ 
lization  systems  described. 

Key  words  (for  information  retrieval) 

Surface  Acoustic  Wave  (SAW)  Resonators,  Aging, 
Drive  Level,  Fabrication,  Electromigracion. 

Introduction 

Systems  requiring  long  term  stability  and 
insensitivity  to  undesirable  effects  provide 
strong  motivation  for  investigating  the  common 
mode  rejection  characteristics  of  appropriately 
configured  pairs  of  surface  acoustic  wave  (SAW) 
resonators. 

Low  absolute  aging  of  Individual  SAW  reso¬ 
nators  depends  on  device  material  composition, 
process  controls  and  cleanliness,  packaging 
environment  and  hermetlclty.  Extended  pre-aging 
appears  to  be  necessary  to  reduce  long  term  drift 
below  1  ppm/yr  with  typical  aging  studies  reported 
in  the  literature  now  approaching  the  two  year 
mark.1_3  An  objective  of  this  study  was  to  iso¬ 
late  as  much  as  possible  those  factors  contribut¬ 
ing  to  the  long  term  instability.  By  a  study  of 
resonator  pairs  in  a  common  enclosure,  effects 
such  as  packaging  could  be  eliminated  or  deduced 
by  evaluating  differential  aging. 


This  paper  will  examine  the  absolute  and  dif¬ 
ferential  aging  of  SAW  resonators  fabricated  as 
pairs  and  sealed  under  vacuum  in  a  common  TO-5 
enclosure.  The  design  of  the  fabrication  process 
and  control  of  the  process  is  described.  Also, 
the  effects  of  power  level  at  which  the  devices 
are  driven  and  the  damage  encountered  under  sus¬ 
tained  conditions  of  high  power  and  elevated 
temperature  are  discussed. 

Fabrication  Process  Design  and  Control 

Reliable  device  performance  requires  a  fabri¬ 
cation  process  that  is  reproducible.  Analysis 
procedures  were  instituted  at  appropriate  steps  in 
the  process  to  monitor  the  condition  of  the  SAW 
substrate  as  it  proceeded  from  wafer  to  chip  and 
finally  to  packaged  environment.  Figure  1  is  a 
flow  chart  of  the  fabrication  process.  As  pilot 
wafers  progress  through  the  fabrication  process 
they  are  subjected  to  a  variety  of  surface  and 
bulk  analytical  measurements  to  determine  the 
presence  of  contaminants,  interdiffusion  of 
deposited  metals,  or  surface  chemical  reactions 
that  may  have  significant  consequences  for  device 
long  terra  stability  after  the  package  is  sealed. * 
This  required  vacuum  baking  prior  to  sealing,  to 
pump  off  moisture  and  accelerate  short  term  sur¬ 
face  reactions.  .Also  a  small  plasma  etch  chamber 
for  frequency  calibration  was  installed  directly 
in  the  dry  nitrogen  glove  box  containing  the 
resistance  weld  sealer.  After  vacuum  sealing  the 
packages  were  bombed  with  helium  and  leak  tested. 
Typically  better  than  90?  of  all  sealed  devices 
indicated  leak  rates  of  less  than  1  *  10“15  std 
cc/sec  air  equivalent. 

Pre-Aging  System 

An  oven  was  designed  for  monitoring  the  aging 
of  36  SAW  resonators  simultaneously.  By  using 
Interchangeable  PC  boards,  either  18  TO-5  packages 
with  resonator  pairs  or  36  packages  containing 
single  resonators  could  be  monitored.  Resonator 
gtoups  of  six  were  connected  to  a  single  pole  six 
throw  (SP6T)  coaxial  switch.  The  six  switches 
were  connected  to  a  seventh  in  a  switching  matrix 
arrangement.  Figure  2  shows  photograph  of  the 
oven,  containing  the  PC  board  into  which  the  T0-5 
packages  were  inserted  and  the  coaxial  switches. 


The  resonators  were  periodically  Interrogated 

_ _  through  the  switches  by  a  network  analyzer,  fre- 

•Worx  supported  by  Defense  Advanced  Research  quency  synthesizer  and  counter,  all  of  which  were 

Projects  Agency,  Contract  No.  MDA  903-81-C-0081 .  controlled  by  a  desktop  computer.  The  resonant 
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Figure  2  Photograph  of  oven  showing  SP6T  coaxial 
switches  and  PC  board  for  TO-5  crystal 
enclosures . 


frequencies  were  sequentially  measured,  the  oven 
temperature  recorded,  and  the  time  logged  from  the 
computer’s  real  time  clock.  In  this  manner,  the 
f requency/tecperature  drift  was  tracked  with 
time.  Software  computed  the  differential  drift 
between  packaged  pairs  of  SAW  resonators. 

Absolute  and  Differential  Pre-Aging  Results 

Pre-aging  studies  were  conducted  on  two 
paired  resonator  configurations  and  two  types  of 
electrode  metal izat ions : 

1.  Pairs:  two  separate  chips  in  one  pack¬ 
age;  evaporated  Cr t A1  electrodes. 

2.  Dual  resonators:  two  resonators  on  a 
common  chip;  evaporated  Cr/Al  electrodes. 

3.  Dual  resonators:  same  as  fl2  but  with 
evaporated  Al  +  4.52  Cu  alloy  electrodes. 

Figures  3,  4  and  5  show  the  absolute  and 
differential  aging  for  each  of  the  three  cases 
indicated.  In  all  cases,  frequency  drift  follows 
a  logarithmic  time  dependence  characteristic  of 
chemi sorptive  reactions.  Small  oscillations  in 
the  aging  were  primarily  due  to  fluctuations  in 
oven  temperature  that  occurred  with  a  24  hour 
cycle.  These  frequency  variations  result  from  an 
offset  between  the  oven  operating  temperature 
(70BO  and  the  turning  point  of  the  crystals 
(85®C).  In  the  case  of  paired  Cr/Al  crystals  the 
absolute  aging  rates  typically  showed  between  5 
and  20  ppm  frequency  increases  within  the  first 
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Figure  5  Absolute  and  differential  aging  for  Cu  doped  A1  dual  resonators. 


30  day*,  typical  drift  r*t«*  fell  below 
O.S  ppa/day  within  48  hour*  and  0. 1-0.5  ppm/wk 
•ft*r  three  week*.  Typical  differential  drift 
rate*  fell  below  0.25  ppa/day  within  48  hour*  and 
typically  0. 1-0.2  ppa/wk  after  three  week* 

(5-10  ppa/y*ar).  Fabricating  the  device*  a*  dual 
resonators  on  the  aaae  cryttal.  Fig.  4,  Improved 
the  differential  aging  rates  but  absolute  aging 
rcaalned  high.  Replacing  the  Cr/Al  metallization 
with  copper  alloyed  aluaimn  remarkably  Improved 
the  absolute  as  well  as  differential  aging,  Fig.  5. 
The  total  pre-aging  amounted  to  less  than  2  ppm  and 
drift  rates  less  than  0.1  ppm/mo  were  obtained  after 
only  3  weeks  pre-aging.  As  expected  the  differen¬ 
tial  aging  showed  a  similar  Improvement.  A  more 
accurate  assessment  of  these  rates  requires  either 
individual  oscillator  measurements  or  longer  pre- 
aging  times.  The  former  Is  currently  under 
investigation  using  a  dual  oscillator  circuit. 

Drive  Level  Studies 


The  effect  of  rf  power  delivered  to  the 
resonators  from  a  50  ohm  source  at  the  resonant 
frequency  of  each  device  was  Investigated  and 
compared  for  three  different  types  of  electrodes: 

(1)  evaporated  chrome,  504,  and  aluminum,  10001, 

(2)  rf  sputtered  aluminum  and  (3)  evaporated  4.5% 
copper-aluminum  alloy,  10001.  It  has  been 
reported  that  excessive  power  levels  In  SAW 
resonators  result  in  drive  level  dependent  fre¬ 
quency  shifting  and  degraded  performance  with 

t lme . 5 , 6 

Figure  6  shows  the  frequency  dependence  on 
drive  level  for  eight  resonators  from  each  of  the 
mecallizatlon  groups  described  above.  From  the 
results  of  the  three  types  of  metallized  devices 
studied  It  does  not  appear  that  one  type  provides 
better  frequency  stability  under  varying  drive 
levels  than  any  other.  This  may  Indicate  that  the 
dependence  of  the  device  frequency  Is  controlled 
primarily  by  thermal  heating  due  to  power  dissi¬ 
pation  In  the  bulk  of  the  quartz,  and  not  by  any 
property  of  the  metallization. 

Degradation  at  High  Drive  Levels 

The  equivalent  circuit  series  resonant  resis¬ 
tance  was  typically  45-75  ohms  for  the  375  MHz  SAW 
resonators  studied  In  this  work  and  a  50  ohm  reso¬ 
nator  will  dissipate  all  the  power  delivered  to  It 
from  a  50  ohm  frequency  source.  At  +10  dBm  drive 
level,  In  excess  of  1  Kw/cc  Is  being  dissipated  in 
the  skin  depth  of  the  central  region  of  the  elec¬ 
trode  pattern.  In  view  of  the  fact  that  the  Q  of 
these  devices  ranged  from  12,000  to  20,000,  and 
che  energy  stored  In  the  device  Is  proportional  to 
the  product  of  the  Q  and  the  energy  dissipated, 
large  stored  energy  densities  are  achieved.  It 
comes  as  no  surprise,  then,  that  electromigration 
effects  due  to  large  piezoelectrically  induced 
potentials,  severe  mechanical  distortions  and 
thermal  stresses  are  likely  to  have  a  degenerative 
effect  on  the  structure  of  the  device,  and  partic¬ 
ularly  the  metal  electrodes. 

SAW  resonators  were  driven  for  periods  ex¬ 
ceeding  15  hours  at  drive  levels  of  +10  dBm. 
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Figure  6  Frequency  dependence  of  drive  level  for 
Cr/Al,  Cu-doped  aluminum,  and  sputtered 
aluminum  metallized  resonators. 


Device  degradation  was  found  to  be  closely  related 
to  the  temperature  at  which  the  drive  level  mea¬ 
surements  were  made.  At  room  temperature  little 
or  no  change  in  che  device  characteristics 
occurred.  At  75"C,  however,  significant  damage 
was  observed  after  prolonged  active  aging  at 
+10  dBm.  Figure  7  shows  the  return  loss  magnitude 
and  polar  response  for  a  SAW  resonator  fabricated 
with  a  Cr/Al  electrode  pattern  before  and  after 
15  hours  of  active  aging  at  75“C.  The  appearance 
of  spurious  resonances,  a  shift  In  the  resonant 
frequency  and  an  Increase  In  the  series  resonant 
resistance  indicates  damage  has  occurred  either  in 
the  acoustic  cavity,  the  lnterdlglcal  transducer, 
or  both. 

Scanning  electron  micrographs  of  the  trans¬ 
ducer  region  of  the  SAW  resonator,  taken  afterward 
confirmed  that  damage  had  occurred.  Small 
hillocks,  typical  of  electromigration  damage,  were 
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Figure  7  Return  loss  degradation  for  +10  dBm 
drive  level  at  75cC, 


found.6  In  addition,  "drag"  type  effects  were 
noticeable  at  the  edges  of  the  electrodes,  as  If 
aluminum  had  been  physically  brushed  or  scraped  In 
the  direction  of  acoustic  motion.  Neither  of 
these  artifacts  were  present  In  undamaged  devices. 

Conclusions 

The  absolute  and  differential  aging  of  SAW 
resonators  was  Investigated.  A  fabrication 


process  which  Included  techniques  for  monitoring 
process  conditions  resulted  in  devices  with  low 
absolute  end  differential  drift  rates  after  a 
short  period  of  pre-aging.  Absolute  aging  rates 
less  than  0.1  ppm/oo  after  only  3  weeks  pre-aging 
at  70*C  were  demonstrated.  Low  differential  aging 
rates  were  obtained  using  dual  resonator  crys¬ 
tals.  Drive  level  studies  on  devices  fabricated 
with  different  types  of  aluminum  metallizations 
have  not  conclusively  shown  that  one  type  of 
metallization  Is  superior  to  another  based  upon 
drive  level.  For  the  purpose  of  long  term 
frequency  stability  however,  copper  alloyed 
aluminum  appears  to  be  superior  to  chrome /aluminum 
metals.  Damage  due  to  elevated  drive  levels  for 
extended  periods  was  found  to  be  more  sensitive  to 
operating  temperature  than  the  drive  level  up  to 
+10  dBm. 
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